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ABSTRACT 

We present a near-infrared spectroscopic study of HD 114762B, the latest- type metal-poor compan- 
ion discovered to date and the only ultracool subdwarf with a known metallicity, inferred from the 
primary star to be [Fe/H] = -0.7. We obtained a medium-resolution (R 3800) Keck/OSIRIS 1.18- 
1.40 /im spectrum and a low-resolution (R - 150) IRTF/SpeX 0.8-2.4 /^m spectrum of HD 114762B 
to test atmospheric and evolutionary models for the first time in this mass-metallicity regime. HD 
114762B exhibits spectral features common to both late- type dwarfs and subdwarfs, and we assign it 
a spectral type of d/sdM9 ±1. We use a Monte Carlo technique to fit PHOENIX/GA/^ synthetic 
spectra to the observations, accounting for the coarsely-gridded nature of the models. Fits to the 
entire OSIRIS J-band and to the metal-sensitive J-band atomic absorption features (Fe I, K I, and 
Al I lines) yield model parameters that are most consistent with the metallicity of the primary star 
and the high surface gravity expected of old late-type objects. The effective temperatures and radii 
inferred from the model atmosphere fitting broadly agree with those predicted by the evolutionary 
models of Chabrier & Baraffc, and the model color-absolute magnitude relations accurately predict 
the metallicity of HD 114762B. We conclude that current low-mass, mildly metal-poor atmospheric 
and evolutionary models are mutually consistent for spectral fits to medium-resolution J-band spectra 
of HD 114762B, but are inconsistent for fits to low-resolution near-infrared spectra of mild subdwarfs. 
Finally, we develop a technique for estimating distances to ultracool subdwarfs based on a single 
near-infrared spectrum. We show that this "spectroscopic parallax" method enables distance esti- 
mates accurate to < 10% of parallactic distances for ultracool subdwarfs near the hydrogen burning 
minimum mass. 

Subject headings: stars: fundamental parameters — stars: individual (HD 114762B) — stars: low- 
mass, brown dwarfs — subdwarfs 



1. INTRODUCTION 

A complete understanding of low-mass stars and brown 
dwarfs is an important goal of astrophysics. Low-mass 
stars (0.08 Mq < M < 0.5 M(7,) a re the most numerous 
objects in our galaxy ijLadal [20061 ) and make up nearly 
half of its total mass. Brown dwarfs are objects that form 
like stars but have insufficient mass to sustain core nu- 
clear reactions (M < 0.08 Mq) and represent the link 
between low- mass stars and extra-solar giant planets. 
Objects with spectral types of M7 and later are known 
as "ultracool" objects and encompass the lowest-mass 
stars and brown dwarfs. Our understanding of these 
objects has dramatically improved over the last decade, 
but there remain several fundamental areas of parame- 
ter space that are relatively unexplored, one of which is 
metallicity. 
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Sub-solar metallicity low-mass stars and brown dwarfs 
are referred to as ultracool subdwarfs. Like their solar- 
metallicity counterparts, ultracool subdwarfs have spec- 
tral energy distributions that peak in the near-infrared 
and are dominated by overlapping molecular absorp- 
tion bands. Their lower metallicities, however, dramati- 
cally alter their spectral properties compared to solar- 
metallicity objects and can ultimately provide insight 
into the influence of metallicity on the chemistry, con- 
densate cloud formation, and temperature/pressure pro- 
files in the atmospheres of ultracool objects. Prominent 
spectral changes resulting from a reduced metallicity in- 
clude suppressed H- and /\-band fiuxes (i.e. bluer near- 
infrared colors) due to in creased collis i on-induced ab- 
sorpt ion by H2 (CIA H2; iLinskvl [l969l : iBorvsow et"an 
[1993), a suppression of metal-oxide bands, and an en- 
hancement of metal hydride bands ijBesselll [1982^ . Ul- 
tracool subdwarfs exhibit large space motions consis- 
tent with thick disk or halo kinematics (jBurgasser et al.l 
|2007[ ). indica ting that th ese objects are probably quite 
old (see, e.g.. iHelmil [2001) . 

The number of spectroscopically-confirmed ultracool 
subdwarfs has rapidly increased from the first identifica- 
tion over a decade a go (LH S 377 , with a spectral type of 
sdM7: iMonet et al.]ll99l : [Gi^[l997l ) to the ~ 45 cur- 
rently known, the vast majority of which have been found 
in the past five year s (see Table 7 of iBurgasser et al.l 
l2007t[Burgassei|[2008D . Many discoveries have been made 
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TABLE 1 

Age and Metallicity of HD 114762A: 
Literature Search (1995 - 2008) 



Reference 



Age (Gyr) 



[Fe/H] 



Haywood (2008) 
Holmberg et al. (2007) 
Gonzalez fc h awslWoi) 
Zhane fc Zhao (2006) 
Reddv et ar72006) 
Saffe et aLT20051 
Valenti & Fischer (2005) 
Santos et al. f2004) ^ 
Gratton et al. (|5o3) 
Laws et aTWoof) 
Heit er & LucTHoOS) 
^dakauc ct al. (2002) 
Fulbright (2000) 
^achaume_et_aL (1999) 
Thevenin fc Idiart7l999') 
Clementini et al. (1999) 



Gonzalez (1998) 
Fuhrrnann (T§98) 
Ng fc B crtclli (1998) 
Henry et al. f 19977 



12.4 (2.8) 
10.6 (2.2) 

14.1 
11.3 (3.0) 
11.8 (3.6) 
7.7 (2.3) 



14 (2) 



11 (1) 



14 (2) 

10.8 (0.8) 
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-0.76 
-0.657 (0.051) 
-0.80 (0.10) 
-0.71 

-0.65 (0.02) 
-0.70 (0.04) 

-0.79Fe I 

-0.78 (0.06) 
-0.74 (0.07) 
-0.7 (0.08) 

-0.72lte; -0.56nlte 
-0.66 (0.07) 
-0.6 (0.06) 
-0.71 (-0.05) 



Adopted Values'* 



11 (3) 



-0.71 (0.07) 



^ The adopted age and metallicity represent the mean and 
rms value of the listed quantities. 

through searches of large near-infrared and proper mo- 
tion surveys inclu ding the Deep Near Infrared Survey of 
the Southern Sky (lEpchtein et al.A997). the Tw o Micron 
AU Sky Survey f2MA SS. ISkrutskie et al.ll200l . and the 
Digi tized Sky Survey (iLepine fc Sharall2005[ ). For exam- 
ple, ILepine fc Schold (|2008[ ) recentlv identified 23 new 
M-type ultracool subdwarfs through template matching 
to Sloan Digital Sky Survey optical spectra. Many more 
field ultracool subdwarfs are expected to be revealed with 
the next generation of deep, multi-epoch all-sky surveys 
such as the Pa noramic Survey Te lescope & Rapid Re- 
sponse System ([Kaiser c t al."2002') and the Large Syn- 
optic Survey Telescope (Tyson 2002). 

The absolute mctallicities of ultracool subdwarfs are 
currently unknown. A commonly-used technique for es- 
timating the metalhcities of late-type objects is to fit syn- 
thetic spectra to their red-optical spec tra and to adopt 
the metallicity of the best-fitting model (iSchweitzer et al.l 
[l999t Ophle et al.l[20p : iBurgasser et aLll2007f ). Evolu- 
tionary models have also been used to estimate the mctal- 
licities of ultracool su bdwarfs in co l or-colo r and color- 
magnitude diagrams (IScholz et ahl l2004bl; IBurgasser] 
l200a IDahn et al.l [20081: ISchilbach et al.ll2009h . It is un- 
clear how reliable these methods are, however, because 
models in this low-temperature, low-metallicity regime 
have not been tested. There is thus a growing need for 
the discovery of ultracool subdwarfs with independently 
derived metalhcities. 

Benchmark objects whose fundamental parameters can 
be independently determined enable valuable tests of at- 
mospheric and evolutionary models. Several solar metal- 
licity benchmark systems have already been used for 
such purposes. On e of the best exam ples is the HD 
130948ABC system (iPotter et al.ll2002D . which contains 
two L4 brown dwarfs in a hierarchical triple system with 
a G2V primary. In this case the age and metallicity of 
the primary star along with the luminosities and dynam- 



ical masses of the brown dwarfs have been directly mea- 
sured and used to t est brown dwarf evolutionary models 
(|Dupuv et al.ll20(39l ). At sub-solar metalhcities. however, 
there have not been any benchmark studies using very 
low-mass objects. One possibility would be to use the 
latest-type objects in globular clusters, which comprise 
coeval stellar populations with uniform chemical compo- 
sitions. Their great distances, however, currently prevent 
spectroscopic studies of the lowest-mass members for 
which even photometric detections have proven di fficult 
to ob tain (Richer et al. 2002; Richer 2006; Riche r et all 
|2008| ). Another approach would be to use nearby ultra- 
cool subdwarfs in binary systems as sub-solar metallic- 
ity benchmarks. Companions to field subdwarfs could 
also act as metallicity calibrators for the field ultracool 
subdwarf population. This technique relies on the con- 
servative assumption that the components of binary sys- 
tems have the same metallicity. Ob servations of bina- 
ries in the field (Desidera et al.ll200^ and in the Hyades 
(jPaulson et al.ll2003l ) have revealed that the components 
had differential [Fe/H] abundances of < 0.02-0.04 dex, 
implying that the assumption of identical metallicity is 
indeed reasonable. Unfortunately, most of the known ul- 
tracool subdwarfs are isolated objects and searches for 
low-mass secondaries to metal-poor s tars have yielded 
surprisingly few bona fide companions (lOsorio fc Martini 
l200llChaname fc Gouldll200llRiaz et al.ll2008D . 

HD 114762B is currently the latest-type object known 
that is gravitationally bound to a metal-poor star (sdF9; 
HD 114762 A). It was discovered and confirmed as a 
proper motion companion in a Lick Observatory adap- 
tive optics search for ste har and subste llar companions 
to exoplanet host stars (jPatience et al.l 12002 !) and was 
later redisc overed in a similar study using CFHT AO 
ijChauvin et al. 2 0()J). HD 114762B is located 3'.'3 (128 
AU at 38.7 pc; Ivan LeeuwenI |2007| ) from its primary, 
which also harbors a high-mass planet, brown dwarf, 
or low-mass star, depending on t he system's inclina- 
tion (Mpsini = 11.6 8 ± 0.96 7\/j„„: 'Latham et al.! ll989l : 
iCochran et "al]ll991l : iHai e 1995; Butler et al. 2006). The 
presence of a massive exoplanet candidate has led to 
numerous studies of the primary's metallicity ([Fe/H]) 
and age, which have been well constrained to be -0.71 
± 0.07 dex and 11 ± 3 Gyr, respectively* (Table [I]). 
Age estimates in the literature for HD 114762A are 
based on theoretical isochrones and chromospheric ac- 
tivity levels. This system may also harbor a debris disk 
based on optical polar imetry and ISO far-IR photometry 
(|Saffe fc G6mezll200l . 

A Keck adaptive optics/NIRS PEC J-band sp ectrum 
of HD114762B was presented bv iPatience et all ()2002f) . 
They found inconsistent spectral features that appeared 
similar to a late-M/early L spectrum based on the 
strength of the 1.4 /im steam band, but also resembled 
a mid-M spectrum based on the strengths of the 1.189 
^m Fe I and 1.313 fxni Al I neutral metal lines and on a 
comparison to J-band spectra of other M dwarfs. How- 
ever, the dearth of ultracool subdwarfs known at the 
time of discovery prevented a comparative analysis with 
other metal-poor objects, and the nature of their ob- 
servations made the analysis difficult because slit-based 

* The quoted errors are rms yaluos based on a literature search 
rather than formal uncertainties. 
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spectroscopy with adaptive optics do e s not preserve the 
continuum shape (iGoto et al.l l2002t iGoto et al.l 120031 : 
iMcElwain et al.ll2007[ ). 

Here we present a near-infrared spectroscopic study 
of HD 114762B. Our goals arc threefold: (1) to better 
characterize this object through a comparison to known 
subdwarfs, (2) to test atmospheric models by compar- 
ing the results of the best-fitting models to the known 
metallicity of this object, and (3) to study the consis- 
tency of the predictions from atmospheric and evolu- 
tionary models. For our tests we use th e GAI A grid 
of synthetic spectra (jBrott fc Hauschildtl 120051 . based 
on the PHOENIX atmos pheric model code) a n d the 
IChabrier fc Baraffd l[l997l CB 9 7). IBaraffe et all (fl997l . 
BCAH97), and |BaraffeeLalJ ([19981, BCAH98) low- 
metallicity evolutionary models (all based on the same 
input physics). These models are widely used through- 
out the literature and arc therefore important to validate. 

In §2 we present our observations and the methods for 
extracting the spectra from the reduced data. In §3 we 
describe the results of the fitting of the models to the 
data and the way in which we derive the fundamental 
parameters from the atmospheric and evolutionary mod- 
els. Finally, we summarize our results and discuss the 
implications of our work in §4. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. OSIRIS Spectrum 

We obtained medium- resolution (R ~ 3800; R = 
A/AA) J-band spectra of HD 114762B with the 
OH-Supp ressing InfraRed Integral-field Spectrograph 
fOSIRIS: lLarkin et al.l[200l) using the natural guide star 
adaptive optics system at Keck II on the night of 2008 
Jan 14 UT. The weather was poor with patchy cirrus 
clouds. Using HD 114762A as the natural guide star, we 
obtained four exposures of HD 114762B with an integra- 
tion time of 600 s and one with an integration time of 120 
s for a total integration time of 42 min. Immediately af- 
terwards we observed the AOV star HD 116960 for telluric 
correction with four exposures of 60 s each. We used the 
50 mas pixel scale for these observations. OSIRIS uses a 
lenslet array to sample a rectangular region of sky, pro- 
ducing 16 X 64 individual spectra which are later merged 
into a data cube covering both spatial and dispersion di- 
rections. The data were reduced and wavelength cali- 
brated using the OSIRIS data pipeline (Krabbc 2004). 
We extracted individual spectra from each observation 
by performing a robust third-order polynomial fit to the 
centroids of each target in the dispersion direction for 
both the the x- and y-coordinates. We then performed 
aperture photometry centered on the polynomial fit po- 
sitions at each wavelength slice using an aperture radius 
of 3 pixels (0'.'15) and an inner and outer sky annulus of 
7 pix (0'.'35) and 10 pix (0'.'50), respectively. The typical 
FWHM was 65 mas. 

Individual spectra were combined by first multiplica- 
tively scaling them to the median-combined spectrum 
and then computing the median of the scaled spectra. 
Measurement errors were determined by computing stan- 
dard errors of the median of the scaled spectra. The 
final spectrum was corrected for telluric absorption us- 
ing a g eneral version of the Spextool reduction package 
(jVacca ct al..,2003; .Gushing et al.,.2004ii) . which was orig- 



inally written for the near-infrared spectrograph SpeX at 
IRTF. 

2.2. SpeX Prism Spectrum 
2.2.1. Observations 

We obtained low-resolution (R ^ 150) 0.8-2.4 /im near- 
infrared spectra of HD 114762B using the near-infrared 
SpeX spectrograph (Rayner et al. 2003) in prism mode 
at IRTF on 2008 Jan 28 UT. Gonditions were not photo- 
metric but the seeing was good with a typical FWHM of 
O'.'S throughout the observations over an airmass range 
of 1.013 to 1.001. We chose a sht length of 15" and a sht 
width of 0'.'5. The close an gular separaticin (3"3 ) and the 
high contrast (AJ ~ 7.6. iPatience et aLll2002[ ) between 
the primary and the companion made the observations 
difficult. With a slit orientation in the radial direction 
from the primary we would have needed to position the 
target near one end of the slit in order to avoid con- 
tamination from the primary, thereby losing a fraction 
of the flux from the companion. We therefore chose a 
slit orientation roughly orthogonal to the line connect- 
ing the primary to the companion in order to capture 
the majority of the flux of HD 114762B near the center 
of the slit. Unfortunately a diffraction spike from the 
primary was situated ~ 3" from the companion so our 
observations include a spectrum of a bright diffraction 
spike along with that of the much fainter science target 
(Figure [H top panel). 

The data were taken by nodding on and off the target 
in the ABBA configuration by about 13", and consist of 
11 cycles of 120 s each for a total integration time of 22 
min. We subtracted the raw data in pairs to remove the 
night sky lines, dark current, and bias offset. The sub- 
tracted pairs were then divided by a fiat-field which was 
created by first scaling to a median-combined fiat and 
then taking the median of the scaled frames. A slight 
residual sky signal was apparent across the entire spec- 
trum. We removed the residual sky at each column by 
fitting a line to the median-combined fiux near both ends 
of the slit in the spatial direction and then subtracting off 
the fitted line. This procedure simply flattened the data 
to zero at both ends of the spatial profile. The effect of 
removing the residual sky was small and the extraction of 
the companion spectrum (see below) was done with this 
flattening and without it; the result on the final spec- 
trum was negligible. We accounted for chip defects in 
the SpeX infrared array detector by using the bad pixel 
mask provided in Spextool. 

2.2.2. Extracting the Spectrum 

Extracting the one-dimensional spectrum from the 
two-dimensional reduced data was not straightforward. 
The main difficulty lay in determining the spatial shape 
of the diffraction spike spectrum whose tail the compan- 
ion was immersed. The peak of the diffraction spike was 
~ 3" from that of the companion in the spatial direc- 
tion. Ideally we would like to know the exact functional 
form of the bright tail so that we could simply subtract 
it off at each column. Practically this meant fitting sim- 
ple functions to the tail to determine the one that best 
approximated the data. 

The spatial profile of the diffraction spike was slightly 
asymmetric so we could not use the opposite tail as a 
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Fig. 1. — Extraction of the HD 114762B spectrum from the SpeX/prism low-resolution data. The upper panel shows a single sky- 
subtracted and fiat-fielded two-dimensional spectrum of the companion HD 114762B and the very bright diffraction spike of the primary 
star HD 114762A. The x-axis is the dispersion direction and the y-axis is the spatial direction (the spatial scale is O'.'IS pix~^). The lower 
panel shows the image after fitting and subtracting the diffraction spike. The approximate locations of the J, H, and K bands are labeled. 



model. We first tested a method of spline interpolation 
across the location of the science target. This technique 
was highly sensitive to both the "tension" parameter of 
the spline as well as the choice of anchor points along the 
tail; the result was a noisy and over-subtracted compan- 
ion spectrum. We tried constant values for the tension 
parameter as well as tension values estimated from the 
best-fitting spline interpolation over the same region on 
the opposite tail. Both methods produced poor results. 
We tested a wide variety of standard functions includ- 
ing a Gaussian, a Moffat, and a number of polynomials 
(from 2'^'^ to 6"^ order) in both normal- and log-space. 
The quality of each fit was judged based on a visual in- 
spection of the fit at each column and of the resulting 
subtracted two-dimensional spectrum. The best result 
was a third-order polynomial fit to the logarithm of the 
intensity values. The functional form in linear flux units 
is 

where / is the intensity at each pixel, y is the pixel coor- 
dinate in the spatial direction, and the values of a, b, c, 
and d are the constants that are fit for at each column. 
In each fit the data were median-combined over 11 pixels 
in the dispersion direction to create a higher signal-to- 
noise spatial profile. Anchor points used in the fit were 
centered around the peak of the companion and were ad- 
justed for a small {^^ 3 pixel) shift caused by a slight rota- 
tion of the data on the array. The adjustment was made 
by fitting Gaussians to the diffraction spike spectrum at 
each column, then fitting a line to the peak of the Gaus- 
sians in the dispersion direction, and then finally offset- 
ting the anchor points from the linear fit. This ensured 
that the anchor points did not move with respect to the 
spectrum of the science target. The fit was then scaled 
to the values of the original (non-median-combined) data 
at the central column by a scale factor that accounted for 
small differences in height of the spatial profile. The scal- 
ing was performed using MPFITEXPR, a robust non-linear 



least squares curve fitting routine in IDL IjMarkwardtl 
I2009f ). Finally, the scaled fit was subtracted from the 
tail of the diffraction spike thereby revealing the spatial 
profile of HD 114762B (Figure [H bottom panel). This 
procedure was applied to every column across the entire 
spectrum for each reduced frame. 

One-dimensional spectra were obtained for each indi- 
vidual frame by summing the flux within a width of 5 
pixels (0'.'75) from either side of the peak of the spatial 
profile of the target spectrum at each column. To ensure 
a smooth aperture radius along the dispersion direction, 
the peak of the spatial profile was determined by first 
fitting a Gaussian to the spatial profile at each column 
and then fitting a line in the dispersion direction to the 
peak of each Gaussian fit. To combine the spectra from 
the individual frames the data were scaled to the median- 
combined spectrum over the H-hand and were median- 
combined again. Telluric correction was performed using 
the Sp extool IRTF data reduction and analys is package 
(v.3.4, iVacca et al.ll2003l : ICushing et aT][200l . We per- 
formed several trial extractions to test the infiuence of 
the aperture radius and anchor points of the fitted func- 
tion on the final coaddcd spectrum. The aperture radii 
were tested with values of 1, 3, 5, 7, and 10 pixels. The 
shape of the spectra changed very little, especially for 
aperture radii > 3 pix. The effect of the anchor point 
positions were also tested by examining the extracted 
spectra and the subtracted two-dimensional spectral im- 
ages of the companion. Large changes of 5 pixels in 
the anchor point positions strongly influenced the overall 
shape of the resulting spectra and introduced consider- 
able noise. We chose anchor points that produced the 
best fits based on a visual inspection of the subtracted 
images and on a comparison of extracted spectra using 
different anchor point regions. 

The contrast between the companion and the diffrac- 
tion spike spectra was especially large sliortward of ~ 
1 /im (Figure [5]). In this region small changes in the 
goodness-of-fit of Equation 1 corresponded to large varia- 
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Fig. 2. — Extraction of the HD 114762B spectrum from the 
SpeX/prism low-resolution data. The spatial scale of the infrared 
array is O'.'IS pix~^. The science target is the small bump seen on 
the steeply rising flux from the primary star's diffraction spike. At 
each column we fit a function (red line, Equation 1) to the spatial 
profile of the diffraction spike using anchor points (demarcated 
by gray dashed lines) that were allowed to shift with the peak of 
the diffraction spectrum across the array. The fit was subtracted 
from the bright tail, revealing the spectrum of the companion 
(plotted as diamonds). The contrast between the science target 
and the diffraction spike was worse in the J band than in the K 
band and may have led to a slight oversubtraction of flux at short 
wavelengths, especially below 1 fim. See § 2.2.2 for details of 
the extraction process and the functional form of the fit. 



tions in the resulting subtracted profile. For that reason, 
the companion spectrum is probably not reliable for A < 
1 fim. While general absorption features are apparent, 
the flux levels are probably not trustworthy. 

2.3. Flux Calibration 

We flux-calibrated the OSIRIS and SpeX spectra us- 
ing published photometry. Differential J, H, and Ks 
magnitudes betwee n HD 114762A a nd HP 114762B are 
hsted in Table 4 of lPatience et all (j2002[ ). Their imag- 
ing observations were obtained using Lick Observatory's 
IRCAL camera with adaptive optics. We used the dif- 
ferential magnitudes combined with 2MASS photometry 
of the primary to compute the IRCAL magnitudes of 
HD 114762B: 



.1 



IRCAL 



AJ 



IRCAL 



2MASS 



+ Jc, 



(2) 



where Jcorr is the correction term used to transform the 
2MASS J-band magnitude to the IRCAL photometric 
system. The same procedure was applied to the H and 
Ks bands. The correction term was derived by calcu- 
lating the differential magnitude of an F9 star in both 
photometric systems. The choice of an F9 star is based 
on the spectral type of HD 114762A (sdF9), which allows 
us to avoid a color correction between the filter systems. 
We performed synthetic photometry on the infrared spec- 
trum of HD 27383 (F9V) which we obtained from the 
IRTF spectral library (Rayner et al., submitted). The 
correction term for the J band is 



-2.5 log 



2.5 log 



J;^^Vcgay2MASS(^)^;^ J 



(3) 



Here TIRCAL j^jmass 

are the transmission curves 
of the J filter in the IRCAL and 2MASS systems, re- 
spectively, while /^^ and f^^^'^ are the flux densities of 
the F9 star and Vega. The 2MASS transmission filters 
include the filter profile together with the atmospheric 
transmission profile, so we multiplied the IRCAL filter 
transmissions with an infr ared telluric profile generated 
using ATRAN ()Lordlll992f) . For Vega fluxes throughout 
this work we used the flux-calibrated spectrum of Vega 
provided in Spextool. The correction terms to convert 
from the 2MASS to the IRCAL systems were < 0.01 
mags for all three bands, which is smaller than the er- 
rors in the 2MASS photometry of HD 114762A and much 
smaller than the errors in the IRCAL differential mag- 
nitudes. The corrections are therefore negligible and we 
consider the 2MASS and IRCAL filters to be the same 
systems in this work. The photometry for HD 114762B is 
listed in Table [21 Errors are derived from the published 
magnitude difference errors and the 2MASS photometric 
errors, added in quadrature. 

We flux calibrated the SpeX data by computing the 
multiplicative constant that scaled the spectrum to the 
flux density level that matched the Xs-band magnitude 
listed in Table [21 The Kg band was chosen because any 
systematic errors caused by the extraction of the spec- 
trum would be minimized in this bandpass (see § 2.2.2). 
The flux calibration scaling factor Cfc was computed in a 
Monte Carlo fashion. An artificial spectrum was gener- 
ated by adding to the data noise drawn from a Gaussian 
distribution with a standard deviation equal to the mea- 
surement error at each pixel. The flux calibration scaling 
factor Cfc was then computed for each artificial spectrum 
in the following manner: 



/A/fg-rKs(A)rfA 

/A/,^^T^s(A)dA 



X 10 



-0.4i<-s,, 



(4) 



where f^p is the Monte Carlo-generated spectrum for 

trial i, and f^'^^'^ and Tks{X) are the same as in Equa- 
tion 3. For each trial i, a new /STs.i magnitude was drawn 
from a Gaussian distribution with a mean value equal to 
the Kg-hand magnitude from Table [21 and a standard de- 
viation equal to the photometric uncertainty. The mean 
flux calibration scaling factor Cfc and its error a-^^ were 

obtained from the distribution of Cfc values for 10^ tri- 
als. The resulting values of Cfc and (7^^. were 1.10 and 
0.10, respectively. The error in the final flux calibration 
level is roughly 9% and is dominated by the iCs-band 
photometric uncertainty. 

The OSIRIS spectrum does not cover the entire 
2MASS J band so we used the SpeX spectrum for flux 
calibration. We smoothed the OSIRIS spectrum to the 
resolution of the SpeX data and then scaled it to the 
same flux level. We assume an identical flux calibration 
uncertainty of 9% for the OSIRIS spectrum based on the 
SpeX spectrum. 

2.4. Synthetic Photometry 
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TABLE 2 

HD 114762B Infrared Photometry (2MASS System) 



Type 








J-H [ctj.h) 


H - Ks icTH-Ks) 


J - Ks (aj^Ks) R-ef 


Aperture Phot. 
Synthetic Phot. 


13.74 (0.10) 
13.99 (0.10) 


13.39 (0.10) 
13.45 (0.10) 


13.01 (0.10) 
12.99 (0.10) 


0.35 (0.14) 
0.54 (< 0.01) 


0.38 (0.14) 
0.46 (< 0.01) 


0.73 (0.15) 1, 2 
1.00 (< 0.01) 3 



References. — f l) IPatience et al.l I I2002I) : f2') ICutri et al.l II2003I) : (3) This- 



E 

CM 




1.35 



Fig. 3. — Line identification in HD 114762B from the medium-resolution OSIRIS spectrum (top) and the low-resolution SpeX/prism 
spectrum (bottom). Major molecular bands are from H2O, FeH, and CO, while prominent atomic features originate from K I, Al I, Na I, 
and Fe I. Measurement errors for both spectra are shown as dotted lines. The portion of the SpeX spectrum that may have been altered 
by the spectral extraction technique is demarcated by the gray shaded region (see §2.2.2). 



We computed synthetic magnitudes and colors from 
our flux-calibrated SpeX spectrum and compared them 
to the published photometry to test whether the ex- 
traction technique altered the slope of the final spec- 
trum. Synthetic magnitudes were derived in a Monte 
Carlo fashion incorporating measurement errors and 
flux-calibration errors. For each trial, an artificial spec- 
trum was generated in the manner described in §2.3 and 
was multiplied by a flux-calibration scaling factor which 
was drawn from a Gaussian distribution with a mean 
value of Cfc and a standard deviation of tr^^. . Synthetic 
J, H, and Ks magnitudes were then computed using 
2MASS transmission profiles. This process was repeated 
10-* times. The mean synthetic magnitudes and the stan- 
dard deviations of the resulting distributions are listed in 
Table [21 Synthetic colors and errors were computed in a 
similar manner and are presented in the same table. 

The synthetic colors and magnitudes from our 



SpeX spectrum disagree with those derived from the 
[Patience et al.l ((2002) differential photometry. The dis- 
agreement appears to progressively worsen at shorter 
wavelengths in the form of a shallower spectral slope in 
the SpeX spectrum. It is unclear whether this difference 
originates in the spectral extraction method used in this 
work or whether it originates with the published pho- 
tometry. A greater offset with decreasing wavelength, 
however, is consistent with an oversubtraction of flux 
in our extraction technique of the companion spectrum 
(§2.2.2). The oversubtraction would be negligible for the 
low contrast Ti'-band region but would be significant for 
the high contrast J-band region. We accounted for this 
difference by dividing the SpeX spectrum into three sec- 
tions and adjusting each section's relative flux to match 
the J-, and Xg-band photometry. The sections were 
chosen to approximate the near-infrared bandpasses: a J 
section (1.1-1.33 /.tm), an H section (1.5-1.8 /im), and a 
K section (2.0-2.35 /xm). We left the if-band spectrum 
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stationary as it was probably least affected by the ex- 
traction process. We refer to the resulting spectrum as 
the "shifted SpeX spectrum" throughout this work. 

3. RESULTS 

3.1. Spectroscopic Characterization of HD 114-762B 

The spectra of HD 114762B reveal that it is a late- 
type object with a spectral energy distribution domi- 
nated by strong overlapping molecular absorption bands, 
most notably the deep 1.4 and 1.9 fim H2O steam 
bands. The strongest atomic absorption lines and the 
major molecular absorption bands are labeled in Fig- 
ure [31 The OSIRIS J-band spectrum shows a pseudo- 
continuum that originates mostly from H2O but that also 
cont ains defining FeH band heads a.t 1.1939 and 1.2389 
um (iPhillips et al.lll987l:lJones et al.lll996l : lCushing et all 
120031 iMcLean et al.l I2007D and a iioticea ble Q-branch 
feature at 1.222 ^m (jCushing et aI1l2003[ ). Prominent 
atomic lines are from Fe I at 1.1886/1.1887 (blended) 
and 1.1976 /im, K I at 1.2436 and 1.2526 fim, and Al I 
at 1.3127 and 1.3154 fim. The SpeX spectrum exhibits a 
rich molecular band structur e which includes the Wing- 
Ford FeH band at 0.99 nm (|Wing fc FordI [19691 ) . other 
strong FeH band heads at 1.239 and 1.194 fim, and the 
2.3 /im CO band. Neutral atomic lines are also visible, 
including a blended K I doublet at 1.177 and 1.245 ^m, 
a blended Al I doublet at 1.314 /im, and a deep Na I 
doublet at 1.14 fim. 

In Figure [H we compare the SpeX spectrum of HD 
114762B to spectra of late- type dwarfs, to so-called "mild 
subdwarfs" which are believed to be only slightly metal- 
poor, and to subdwarfs in the top, middle, and bot- 
tom panels, respectively. The spectra are from the 
SpeX Prism Spectral Library^ and are normalized at 
1.275 /im. HD 114762B closely resembles the M9 ob- 
ject 2 MASS 17072343-0558249 ([McElwain fc Bureassed 
l'2006[) between 1.00-1.35 fim, but has a deeper 1.4 fim 
H2O steam band than any of the normal dwarfs. A much 
better fit is obtained when compared to the mild subd- 
warfs. The general shape of the spectral energy distribu- 
tion matches quite well, except in the K band, where the 
spectrum of HD 114762B appears slightly higher than the 
mild subdwarfs. This may be caused by a slight oversub- 
traction at shorter wavelengths (see § 2.2.2) which, when 
normalized to 1.275 fim, would result in a raised i^-band 
fiux. Finally, a comparison to subdwarf templates shows 
poor matches across the entire wavelength range. 

The best overall fit appear s to be to the d/sdM9 o bject 
2MASS 00412179-^3547133 (iBureasser et al.l[200l . We 
therefore adopt a spectral type of d/sdM9 ± 1 for HD 
114762B. In Figurc[S|we compare 2MASS 0041+35 to the 
shifted SpeX spectrum. The agreement is satisfactory 
over the J and K bands, but the H band is significantly 
lower. This reduced H band flux would only worsen the 
agreement between the shifted spectrum and the spectra 
from Figure[4l we therefore suspect that this lowered flux 
is an artifact of shifting the segments of the SpeX data 
to published photometry which had large uncertainties 
(color errors were ~ 0.15 mags). In addition, we note 
that the d/sdM9 spectral type for HD 114762B is con- 
sistent with the sequence of medium-resolution subdwarf 



spectra shown in Figure [SI based on the strength of the 
1.4 fim steam band and the strengths of the 1.2436 and 
1.2526 fim K I lines. 

Another line of evidence independently bolsters the 
metal-poor nature of HD 114762B. This object has a blue 
J - A's color (0.73 mag) compared to normal dwarfs with 
(optical) spectral types later th an ~ M6, which ha ve J - 
Ks values > 1.0 mags (see, e.g.. iKirkpatrickl [20081 ). Blue 
near-infrared colors may be produced through strong 
CIA H2 caused by a reduced metallicity. It can also 
originate from thin or large grain condensate clouds 
(jBurgasser et al. l2008aD. Clouds begin to appear in th e 



^ Maintained by Adam 

jhttp: / /www.browndwarfs.org/spexprismj 



Burgasser 



at 



latest" M dwarfs ([Jones fc Tsuii|[l997riAllard et al.|[200l[) , 
but their influence on the near-infrared spectral slope 
only becomes significant in L dwarfs, as is evident by the 
increased spread in J - Kg color in that spectral class. 
The blue near-infrared colors observed in HD 114762B 
arc therefore indicative of a reduced metallicity. 

3.2. Physical Parameters From Atmospheric Models 

There is no standard method for fitting a coarsely- 
sampled grid of synthetic spectra (e.g., A[M/H] ~ 0.5, 
ATcff ~ 100 K, and Alog g ~ 0.5) to observations. Gen- 
eral techniques for fitting models to data usually rely on 
minimizing (some form of) the statistic, equal to 
ifi ~ •^i)^/"'!' where fi, Aii, and ai are the data, the 
model prediction, and the measurement error at point i, 
respectively. When the random variable {fi - A4i)/ai is a 
standard normal deviate, the statistic approximates 
a distribution. In that case confldence regions for 
the best-fitting parameters may be approxim ated by cal- 
culat ing the Xniin"Pl^^"'^'^i^^*^iit value (e.g., [Press et all 
[2OOI . 

Atmospheric models have known systematic errors in 
the form of incomplete opacity sources, in which case 
the method of error analysis may not be applicable. 
Thus the question arises of how correct the models must 
be for this technique to remain reliable. Small devia- 
tions from "truth" in the model will have little influence 
on the resulting confldence regions, but large deviations 
will make this method more and more incorrect. One 
way to assess this applicability is to study the normal- 
ity of the distribution of deviates {fi - M.i)/ai for the 
best-fltting model. A near-perfect model will have devi- 
ates that are approximately normally distributed, but an 
inadequate model will have non-normally distributed de- 
viates. Note that the condition of normally distributed 
deviates does not guarantee that the model is good, but 
a good model should have normally distributed deviates. 
A test for normality of the deviates may therefore serve 
as a useful assessment of the correctness of a best-fltting 
model. Reduced x^ values are also used as a goodness-of- 
flt measure, but they rely on an accurate determination 
of measurement errors, a task that is difflcult to correctly 
perform. 

3.2.1. Method 

In this study we first determine the synthetic spec- 
trum that best fits the data using a Monte Carlo ap- 
proach described below. We then make use of two pa- 
rameters to determine the goodness-of-fit: the reduced 
X^ value and the signi ficance level of the D'A gostino- 
Pearson normality test ([D'Agostino et al.|[1990f) applied 
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Fig. 4. — Comparison of HD 114762B (black) to spectra of latc-typc objects using the same instrument at IRTF (SpeX/prism). The top 
panel is a comparison to late-type dwari^s, the middle to mild subdwarfs, and the bottom to subdwarfs. The deep 1.4 and 1.9 /im H2O 
steam bands are poorly matched by the dwarfs and the subdwarfs. We assign a spectral type d/sdM9 it 1 to HD 114762B based on its close 
resemblance to the mildly m etal-poor object 2MAS S 0041-(-35. The d ata ar e from the SpeX Prism Spectr al Library and the observations 
were originally presented bv lBurgasser et al.l l|2004l ). [Burgasserl l|2004l ) . and IMcElwain fc Burgasserl l|2006l ^. The full names of the objects 
used for comparison are, from top to bottom panels, 2MASS 04071296+1710474, 2MASS 00552554+4130184, 2MASS 17072343-0558249, 
2MASS 13593574+3031039, 2MASS 15590462-0356280, 2MASS00412179+3547133, LHS 377, 2MASS 01423153+0523285, SSSPM 
1013—1356, and 2MASS 16262034+3925190. The portion of the SpeX spectrum that may have been altered by the spectral extraction 
technique is demarcated by the gray shaded region (see §2.2.2). 




1.0 1.2 1.4 1.5 1.8 2.0 2.2 



Fig. 5. — S peX spectrum shifted to the photometry of 
IPatience et al.l |2002l . solid black) and compared to 2MASS 
0041+35 (d/sdM9, gray), which is the best-fitting reference 
object to the original unshifted SpeX spectrum (dotted black). 
All spectra are normalized to 1.275 fiia. The i^f-band flux of 
the shifted SpeX spectrum is significantly lower than that of 
2MASS 0041+35, which is probably a result of scaling the spectral 
segments of HD 114762B to the published photometry which had 
large errors. 



to the best-fitting model deviates. The latter test com- 
pares the skewness and kurtosis of a distribution to that 
of a normal distribution. If the deviates are normally 
distributed and the reduced value is low (we impose a 
cutoff of 2.0) then the model is deemed "good" and ap- 
proximately Gaussian confidence limits can be derived 
using the aforementioned Xmin"Pl'^S"'^0'^st=''^t method. If 
these criteria are not satisfied then we simply report the 
physical parameters of the top few best-fitting models 
without errors. 

Wc use the GAIA grid of synthetic spectra, derived 
from the P HOENIX stellar atmospher e code (GAIA ver- 
sion 2.6.1, iBrott fc Hauschildd [20051 : iHauschildt et all 
|1999[ ). as models for our observations. We made use of 
4030 models with zero a-element enrichment compared 
to solar values covering the cool stellar and substellar 
regime: -4.0 < [M/H] < 0.5 (A[M/H] = 0.5), 2000 K 
< Teff < 5000 K (ATeff = 100 K), -0.5 < log g < 5.5 
(Alog g = 0.5). Examples of these models over a range 
of metallicities, temperatures, and surface gravities are 
presented in Figure [T] 

Our approach for fitting the synthetic spectra to the 
observations makes use of a Monte Carlo method previ- 
ously used by sever al authors (e.g.. iSaumon et "al]l2006l . 
ICnshing et atllM)! . The GAIA models were Gaussian 
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LHS 3409 sdM4.5 




1.15 1.20 1.25 1.30 1.35 

\ (/im) 

Fig. 6. — J-band spectra of increasingly later- type metal-poor 
objects compared to our OSIRIS spectrum of HD 114762B. 
The data for LHS 3409, LH S 1135, and LSR 2036-1-5059 are 
from ICushing fc Vaccal I I2006I) taken with SpeX in SXD mode 
at IRTF. Spectra are normalized to 1.27 fira and are shifted 
by a constant. HD 114762B is consistent with a spectral type 
later than LSR 2036-f5059 (sdM7.5) based on the depth of 
the 1.2436 and 1.2526 fim K I lines and the strength of the 
1.4 /im H2O steam band. The spectrum of HD 114762B has 
been Gaussian smoothed to the same resolution as the other data. 



smoothed to the resolution of the SpeX and OSIRIS data 
and were resampled onto the same wavelength grids. The 
synthetic spectra provide the emergent flux density at the 
surface of a star and have to be scaled to the flux density 
observed at Earth. For each model k we compute the 
scaling factor Ck by minimizing the statistic 



X 



E 



fi — CkSkA 



(5) 



where fi and ct^ are the observed flux density and the 
measurement error, respectively, £]^_i is the emergent flux 
density for model k and pixel i, and n is the number of 
data points used in the fit. The scaling factor is given by 



Ef^£k,^/cTf 



(6) 



and is related to the stellar radius R and the distance d 
to the source by Cfc = (R/d)^. 

The GAIA model with the lowest value may change 
for different observations of the same source as a result of 
random noise in the data. To account for this uncertainty 
we generated a large set of simulated data and then de- 
termined the best-fitting model for each new spectrum. 
The simulated data were created by adding noise to the 
original spectrum using the measurement errors. To ac- 
count for the uncertainty in the fiux calibration, which is 
important for the error in the model scaling factor a-^^ , 
the artificial spectrum was adjusted by a flux calibration 
scaling factor drawn from a Gaussian distribution with a 
mean Cfc and a standard deviation a-^^ . Following the 
approach of ICushing et al.l l)2008f ). the 20 models with 
the lowest original values were then refit to the arti- 
ficial data sets for 10,000 trials of simulated spectra. We 
adopted the model with the highest fraction of val- 



ues as the best-fitting model; we re fer to that fraction a s 
/mc, the "Monte Carlo fraction" ijCushing et al.ll2008f) . 
A mean scaling factor Ck was calculated by averaging 
the individual Ck values for trials where the model with 
the highest resulting Monte Carlo fraction had the lowest 
value. The scaling factor uncertainty cr^^ represents 
the standard deviation of the distribution of Ck values. 

3.2.2. Results 

Accurate atmospheric models should yield similar best- 
fitting physical properties for fits to observations of the 
same target at different wavelengths and spectral resolu- 
tions. We fit the GAIA synthetic spectra to the OSIRIS 
spectrum (from 1.18-1.35 /im) and to the SpeX spectrum 
(from 1.00-2.35 ^m; Figure E]). We also fit the shifted 
SpeX spectrum (Figure [H see § 2.4) as well as the in- 
dividual near-infrared bandpasses of the original SpeX 
spectrum (from 1.10-1.33 /im for J, 1.50-1.80 /im for H , 
and from 2.05-2.35 /im for A') to study the wavelength 
dependence of the fits. Our Monte Carlo fitting of the 
shifted SpeX spectrum incorporated measurement errors 
from both the observed spectrum and the published pho- 
tometry. For the shifted SpeX spectrum, each spectral 
segment was adjusted by a scaling factor drawn from a 
normal distribution with a standard deviation equal to 
the error in the published photometry. Then measure- 
ment errors were drawn and added back to the data. 
Results from the fits are presented in Table [3l In Figures 
[TUl and [11] we present contour plots of slices of the 
cubes near the best-fitting models for fits to the entire 
OSIRIS and SpeX spectra. 

The resulting metallicities and surface gravities vary 
greatly among the best-fitting models from Table [H but 
there is good agreement in the derived effective temper- 
atures (2500-2800 K). The fitting results that suggest a 
low surface gravity for this object are inconsistent with 
the known high sur face gravities of old low-mass stars 
an d brown dwarfs (jLopez-MoralesI 120071 : see Figure 2 
of iJao et all 120081 for a graphical version of the data) . 
The fit to the OSIRIS spectrum resulted in a best-fitting 
model with [M/H] = -0.5 and log g = 5.5, in accord 
with the metallicity of the primary^ and the high surface 
gravities of late-type objects. Reduced x^ values are gen- 
erally quite high and the significance levels (p) from the 
D'Agostino-Pearson normality tests are mostly low, in- 
dicating that the best-fitting models do not match the 
observations well. A visual inspection confirms this, al- 
though visually the OSIRIS fit appears to be quite good 
( "chi- by-eye" ) . If the p value is lower than 0.005 then we 
reject the null hypothesis that the deviates are normally 
distributed. In no case were both p > 0.005 and Xmin/^ 
< 2.0 so a formal error analysis on the fitted physical 
parameters is not performed. 

Fits to the SpeX spectrum yielded a large range of 
metallicities (-2.0 < [M/H] < -0.5 dex) and surface 
gravities (1.0 < log g < 5.0). It is unclear whether 
these discrepancies result from a failure of the GAIA 
models or originate from the SpeX data themselves, e.g. 
as a result of the atypical spectral extraction technique 

® We do not different iate between [Fe/H ] and [M /H] for our tests, 
although we note that IValenti &: Fischeil I|2005I 1 derive a value of 
[M/H] = -0.52 for HD 114762A. This difference does not affect the 
analysis nor the results, however. 
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[0.0/2700/5.0 
-1.0/2700/5.0' 
-2.0/2700/5.0 



-0.5/3000/5.0 
-0.5/2700/5.0 
-0.5/2400/5.0 




-0.5/2700/5.5 

-0.5/2700/4.5' 
-0.5/2700/3.5 




Fig. 7. — The PHOENIX/GA/yl v. 2. 6.1 atmospheric models I IBrott Hauschif3tll2005l) showing the variations in emergent flux density 
with metaUicity (top), temperature (middle), and surface gravity (bottom). The result of a reduced mctallicity includes a suppression of 
H and iiT-band flux from collision-induced H2 absorption, producing bluer near-infrared colors. The model notation is as follows: [[M/H], 
Teff (K), log g]. 



^ 0.2 



[-0.5/2800/5.5] 




Fig. 8.— Resuhs of fltting the OSIRIS (top) and SpeX (bot- 
tom) spectra. The best-fltting GAIA model to the 1.18-1.35 fim 
region of the OSIRIS spectrum was [-0.5/2800/5.5]. The best- 
fitting model to the SpeX spectrum from 1.00-2.35 /im was [- 
2.0/2500/1.0]. The metaUicity from fitting the OSIRIS spectrum 
([M/H] = -0.5) is consistent with that of the primary star HD 
114762A ([Fc/H] = -0.70). The gray regions at the bottom in- 
dicate the sections of the observations that were used in the fit. 



(§2.2.2). To address this question we fit tlic models to 
published near-infrared spectra of four ultracool subd- 
warfs with similar spectral types to that of HD 114762B. 
If the models predict mutually consistent values with 




Original SpeX 
Shifted SpeX 



[0.5/2600/3.0] 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Fig. 9. — The best-fitting model (red) to the shifted SpeX spec- 
trum (solid black) compared to the original unshifted SpeX spec- 
trum (dotted black). The resulting metaUicity ([M/H] = -1-0.5) is 
inconsistent with the value of the primary star ([Fe/H] = -0.70). 
The gray regions at the bottom indicate the sections of the obser- 
vations that were used in the fit. 



only small variations in each parameter then that would 
indicate that it is our SpeX spectrum that may be 
skewing the results of the fits. We examined the spec- 
tra of three mildly metal-poor d/sdJVIS objects (21VIASS 
01151621+3130061, 2JVIASS 15561873+1300527, 2JVIASS 
15590462-0356280; NIR spectral types) and one d/sdM9 
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TABLE 3 
Best-Fit PHOENIX/GA/A Models 



Spectral Region 

Used in Fit [M/H] T,a (K) log g /mc P xl,J^ 



OSIRIS 



1.18 - 1.35 Aim 


-0.5 


2800 


5.5 


0.998=' 


< 10-3 


7.29 


Absorption Lines'" 


-0.5 


2600 


5.0 


0.824=1 


< 10-=* 


5.58 




0.0 


2500 


5.0 


0.176 


< 10^=* 


5.69 


SpeX 


1.00 - 2.35 lira 


-2.0 


2500 


1.0 


0.500=" 


< 10-3 


6.47 




-2.5 


2500 


0.5 


0.343 


< 10-3 


6.49 




-2.0 


2600 


1.0 


0.157 


< 10-3 


6.53 


1.00 - 1.35 Atm 


-1.5 


2600 


3.5 


0.697=" 


0.104 


3.82 


-1.0 


2600 


4.0 


0.231 


0.305 


3.89 


1.50 - 1.80 Atm 


-1.5 


2600 


1.5 


0.422=" 


0.911 


5.87 




-1.0 


2600 


2.0 


0.378 


0.921 


5.87 




-3.0 


2700 


-0.5 


0.174 


0.902 


5.99 


2.05 - 2.35 /im 


-0.5 


2500 


5.0 


0.527=" 


< 10-3 


1.47 




-0.5 


2600 


5.0 


0.287 


< 10-3 


1.48 


SpeX Shifted to Match Published Photometry'^ 


J/H/K bands 


0.5 


2600 


3.0 


0.392 


< 10-3 


8.09 




0.5 


2700 


3.5 


0.130 =" 


< 10-3 


7.20 



Note. — Synthetic spectra are from the PHOENIX/GA/A v. 2. 6.1 model atmo- 
sphere code ]Brott 8z Hauschildt|[2005D . Column 1 lists the region of the spectrum 
used in the fit. The resulting metallicities, effective temperatures, and surface grav- 
ities are hsted in Columns 2, 3, and 4, respectively. The Monte Carlo fraction (/mc) 
is listed in Column 5 (§3.2.1). Models with /mc ^ 0.1 are reported, and the model 
with the highest fraction is shown in bold. A D'Agostino-Pearson normality test is 
applied to the deviates for each model. The significance level p of the null hypoth- 
esis, that the deviates are normally distributed, is listed in Column 6. The null 
hypothesis is rejected for values of p below 0.005. Reduced values for fits to the 
observed data rather than the simulated data are listed in Column 6, where i/ — n 
- m — 1 for n data points and m parameters used in the model. Note that m — 4 
for all fits except for fits to the absorption lines, in which case m — 7 (Tcff , [M/H], 
log 5, and one scaling factor for each of the four wavelength ranges). 

^ This model had the lowest value in a fit to the observed data. 

^ Five lines were used: Fe I (1.1886/1.1887 tiui: blended), K I (1.2436 ^m), K I 

(1.2526 Aim), Al I (1.3127 ^m), Al I (1.3 154 ^m). 

Shifted to match photometry from iPatience et al.1 (|20Q2|^ . The fit was 
performed to the J, H, and Kq bands, skipping over the steam bands. See §2.4 for 
details. 



obje ct (2MASS 0041 2 179+ 3547133; NIR spectral type) 
from lBurgasser et al.l ()2004l ). made available through the 
SpeX Prism Spectral Library. We used the measurement 
errors from our SpeX observations scaled to the square 
root of the ratio of the respective integration times to 
determine /mc values^. We first fit the entire 0.80-2.35 
/xm range of the prism spectra, presented in Figure fT2l 
The resulting metallicities were 0.5, 0.0, 0.5, and -3.0 
dex for the four (presumably) mildly metal-poor objects. 
To test the wavelength dependence of the predictions 
we also fit the mild subdwarfs from 1.00-2.35 /.tm and 
over the individual J, if, and K bandpasses (Table |4]). 
The resulting metallicities and surface gravities varied 
considerably among the four objects. Fits to different 
spectral regions of the same object also produced highly 
discrepant results. The effective temperatures were rela- 
tively constant, however, with most fits giving 2500 K or 

Initially we used the errors provided in the SpeX Prism Spec- 
tral Library which were derived from the observations, but reduced 
chi-squared values were quite high, indicating that the errors were 
likely underestimated. 



2600 K. The overall inconsistency implies that the metal- 
licities and surface gravities resulting from fits to low- 
resolution near-infrared spectra of ultracool subdwarfs are 
probably not reliable. 

We also fit the metal-sensitive absorption lines in our 
OSIRIS spectrum to determine whether this fitting tech- 
nique was superior to fitting the entire J band spectral 
region, which emphasizes the overall shape of the spec- 
trum rather than the strengths of the individual absorp- 
tion lines. We identified prominent lines in the OSIRIS 
observations which were sensitive to metallicity, temper- 
ature, and surface gravity in the GAIA models (sec Fig- 
ure [T31). Five lines were chosen: Fe I (1.1887 /im), K I 
(1.2436 /im), K I (1.2526 /im), Al I (1.3127 /im), and 
Al I (1.3154 /.tm). We fit these lines in the OSIRIS data 
using a fitting range of ± 30 A from the center of each ab- 
sorption line, except for the aluminum doublet, for which 
the range was —30 A from the 1.3127 /im line to -f30 A 
from the 1.3154 /xm line. The data were normalized to 
unity at the shortest wavelength of the fitting range for 
each line and a separate scaling factor was computed for 
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TABLE 4 

PHOENIX/GA/yl Model Fits to Published Late-Type Mild Subdwarf SpeX 

Spectra 



Object- [M/H] Teff (K) log g /mc P X^J^ 



Fit from 0.80-2.35 /im 



2MASS 0115-1-31 (d/sdM8) 


0.5 


2600 


4.0 


0.999 


< 10-=* 


6.95 


2MASS 1556-1-13 (d/sdM8) 


0.0 


2500 


3.0 


1.000 


< 10-=* 


9.24 


2MASS 1559-03 (d/sdM8) 


0.5 


2600 


3.5 


0.814 


< 10-3 


7.24 


2MASS 0041-1-35 (d/sdM9) 


-3.0 


2300 


0.0 


0.972 


< 10^3 


15.04 




Fit from 1.00-2.35 


^m 








2MASS 0115-1-31 


0.0 


2700 


4.5 


1.000 


< 10-=* 


3.84 


2MASS 1556-1-13 


-1.5 


2600 


2.0 


0.959 


< 10-3 


4.66 


ziVlAoo iooy— Uo 


0.5 


2600 


A 

4.0 


U.ooy 


^ 1 n — 3 


D.Do 


oA/r A oc nn/11 i qc^ 
ZMAoo UU4i-t-iiO 


-1.5 


2600 


Q n 
o.U 


n non 


^ 1 n — 3 
< iU 


O.oD 




Fit from 1.00-1.35 


/^m 








2MASS 0115-1-31 


-1.5 


2600 


3.0 


0.231 


0.206 


3.47 


2MASS 1556-1-13 


-2.5 


2600 


1.5 


0.398 


0.484 


3.21 


2MASS 1559-03 


-1.5 


2600 


2.5 


0.637 


0.866 


3.28 


2MASS 0041-1-35 


-1.0 


2300 


4.0 


0.373 


0.404 


5.26 




Fit from 1.50-1.80 


^m 








2MASS 0115-1-31 


-0.5 


2700 


4.0 


0.822 


0.060 


3.25 


2MASS 1556-1-13 


-2.0 


2600 


1.0 


0.852 


0.920 


4.68 


2MASS 1559-03 


0.0 


2700 


4.5 


0.429 


0.664 


6.48 


2MASS 0041-1-35 


0.0 


2500 


5.0 


0.782 


0.702 


5.38 




Fit from 2.05-2.35 


^m 








2MASS 0115-1-31 


0.0 


2800 


5.0 


0.498 


0.865 


1.67 


2MASS 1556-1-13 


0.5 


2800 


5.5 


0.855 


0.022 


2.17 


2MASS 1559-03 


-0.5 


2700 


4.0 


0.370 


0.708 


4.24 


2MASS 0041-1-35 


-1.0 


2700 


4.5 


0.487 


< 10-3 


2.56 



" The full names of the objects are 2MASS 011.51621+3130061, 2MASS 
15561873+1300527, 2MASS 15590462-0356280, and 2MASS 00412179+3547133. 




Fig. 10. — Contour plots of values from fitting the GAIA syn- 
thetic spectra to our OSIR,IS spectrum. The top row shows slices in 
metallicity near the global minimum; the bottom row shows slices 
in surface gravity. The global minimum is located at [M/H] = -0.5, 
Tefi = 2800 K, log g = 5.5. The contour levels represent 1.02, 1.3, 
1.9, 2.5, 3.1, and 3.7 times the . value, where . = 8222.7 

' ' ' -^mm ' '^min 

(dof = 1128). Although confidence regions of contours cannot be 
computed using the Xmin-plus-constant method (see § 3.2.1), the 
overall shapes are useful aids for visualizing how constrained each 
parameter is. Local minima can also be easily spotted in this way. 
In this case the temperature is well constrained, but the surface 
gravity and metallicity are not. 




5000 4000 3000 4000 5000 4000 

T„, (K) 

Fig. 11. — Contour plot s of values from fits to our SpeX 
spectrum (similar to Figure FlOt . The global minimum is located at 
[M/H] = -2.0, Tcft = 2500 K, log g = 1.0. The contours represent 
values of 1.02, 2.0, 4.0, 6.0, 8.0, and 10.0 times the . value, 

11/11 ■^•-niiii ' 

where xl,in = 2324.4 (dof = 359). 



each of the four spectral segments. Each model was then 
fit independently to the four regions and the resulting 
values were summed to produce a final value for 
that model. The parameters from the best-fitting mod- 
els using this technique are similar to those for fits to the 
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2UASS 1556+13 (d/sdM8) 
[0.0/2500/3.0] 




Fig. 12. — Spectral fits to SpcX/prism d ata for objec ts with sim- 
ilar spectral types to HD 114762B (from [Burgasscr et al. 2004!). 
The metallicities and surface gravities from the best-fitting mod- 
els vary considerably among the objects for fits to the 0.80- 
2.35 fim region. The full names of the mild subdwarfs are 
2MASS J01151621+3130061, 2MASS J15561873-I-130027, 2MASS 
J15590462-0356280, and 2MASS J00412179-I-3547133. The gray 
lines at the bottom indicate the sections of the observations that 
were used in the fit. 



entire OSIRIS J band spectrum f Figure [M]) . although in 
this case the effective temperature is reduced by 200 K 
to 2600 K^. 

The effective temperatures of the best-fitting GAIA 
models (2500-2800 K) are systematically hotter than 
the temperatures of other late M-type objects which 
have been der i ved u sing other techniques. For example, 
iGautier et al.l (|2007f ) use 24 fim photometry from Spitzer 
Space Telescope to derive effective temperatures of M 
dwarfs using the infrared flux method. They find effec- 
tive temperatures of 2150-2450 K for spectral typ es be- 
tween M8 and M9. 5. Similarlv. iGolimowski et all (|2004f ) 
derive effective temperatures between 2000-2525 for spec- 
tral types between MB. 5 and M9.5 based on LBoi-Ieff 
relations from evolutionary models. The hotter effec- 
tive temperatures we inferred from this study may have 
two (non-mutually exclusive) causes: systematic errors in 
the models or different temperature scales caused by the 
metal-poor nature of HD 114762B. The GAIA models 
do not include the effects of dust, which can be a signif- 
icant source of continuum opacit y in late M-type and 
L-type objects (iTsuii et al." '19961 : iJoiies k Tsuii' 119971: 
Leinert et all 120001 : lAUard ct a!T[200H iPavlenko et all 



2006D . 



° We also tried fitting the models using the pseudo-equivalent 
widths of the model absorption lines relative to those of the data. 
We calculated a chi squared statistic 



2 

Xew 



E 



EW 



M,l 



(7) 



over lines i = 1 to 5. Here EWf^i is the pseudo-equivalent width of 
the observed data for line I, and EWj^^j is the pseudo-equivalent 
width of the model line. A Monte Carlo simulation was used to 
determine Monte Carlo fractions, /mc values were quite low, in- 
dicating that this method did a poor job of distinguishing among 
the best-fitting models. The results were also dependent on the 
way the pseudo-equivalent widths were calculated, as the pseudo- 
continuum is far from smooth. This method seems to be inferior 
to fits over the absorption lines. 



The AMES-Dusty grid (jAllard et al.ll2001f) is a popular 
alternative set of low-temperature models which include 
the limiting effects of dust, but the sub-solar metallic- 
ity AMES-Dusty models only extend as cool as 2800 K. 
Nevertheless, to gauge the differences between non-dusty 
and dusty models, we fit the solar-mctallicity AMES- 
Dusty models (AToff = 100 K; Alog g = 0.5) to our J- 
band OSIRIS and 1.00-2.35 /im SpeX/prism spectra of 
HD 114762B. The best-fitting AMES-dusty model to the 
OSIRIS spectrum is [2800/5.5] and to the SpeX spectrum 
is [2600/4.0]. These effective temperatures arc similar to 
those of the best-fitting GAIA models and so the use 
of non-dusty models does not appear to be the primary 
cause of the hotter temperatures. However, it remains 
possible that both sets of models produce systematically 
hot effective temperatures. 

An alternative explanation for the hotter effective tem- 
peratures is the metal-poor nature of HD 114762B com- 
pared to the solar-metallicity M dwarfs used in previous 
studies. As we discuss in § 3.2.4, ultracool subdwarfs 
appear to be overluminous compared to field ultracool 
dwarfs of the same spectral class, possibly a result of 
a hotter effective temperature scale for ultracool subd- 
warfs. A systematic study of ultracool subdwarf effective 
temperatures using other methods, like the infrared flux 
method, will provide insight into this discrepancy. 

3.2.3. The Luminosity of HD 1147623 

We calculated the luminosity of HD 114762B using its 
bolometric flux and the (revised) HIPPARCHOS paral- 
lax measurem ent of HD 114762A (25.87 ib 0.76 mas, o r 
38.651};}^ pc: lvan Leeuwenll2007l: [Perrvman et al.lll997f) . 
The bolometric flux was computed twice using the origi- 
nal and shifted SpeX spectra to test whether shifting the 
spectrum significantly influenced the resulting luminos- 
ity. To determine the bolometric flux of the unshiftcd 
data we first created an artificial spectrum by adding 
noise drawn from the measurement errors. We then used 
the best-fitting model as a bolometric correction by at- 
taching a short-wavelength segment (10~^ jjjiv < A < 
1.00 ^m) and a long- wavelength segment (2.35 /xm < A 
< lO'^ /im) to the artificial spectrum. To account for 
the uncertainty in the flux calibration we adjusted the 
spectrum by a scaling factor drawn from a Gaussian dis- 
tribution with a mean Cfc and error cr^^ . The bolomet- 
ric flux was then computed by integrating the spectrum. 
This process was repeated 10,000 times, from which a 
mean bolometric flux and error were determined. A sim- 
ilar process was carried out for the shifted spectrum. 
Each artificial shifted spectrum consisted of J, H , and K 
spectral regions to which noise was added. In addition, 
short-wavelength (10^"^ /xm < A < 1.10 /im) and long- 
wavelength (2.35 /im < A < 10'^ /xm) segments from the 
best-fitting model were attached to the artificial data, as 
were model contributions between the J to H and H to 
K filter bandpasses. The spectrum was multiplicatively 
scaled by a Monte Carlo-generated flux-calibration scal- 
ing factor and the bolometric flux was computed. This 
process was repeated 10,000 times. 

The uncertainty in the bolometric luminosity incorpo- 
rates the error in the parallactic distance measurement 
((Td) and the error in the bolometric flux (ctfboi) in the 
following way: 
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Fig. 13. — Relative strengths of the J-band absorption lines in the GAIA model atmospheres for changing physical parameters. The 
depth of the lines are sensitive to changes in metallicity (top), temperature, (middle), and surface gravity (bottom). The combination of 
these lines provides a good estimator of all three parameters when fitting synthetic spectra to observations. The models arc smoothed for 
better rendering. 
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-0.5/2600/5.0] 
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Fig. 14. — Results from fitting the OSIRIS data to the absorp- 
tion lines. The metallicity and surface gravity are consistent with 
the metallicity of the primary star ([Fe/H] = —0.70) and the high 
surface gravity of late-type objects. The effective temperature is 
consistent with that predicted by the CB97 evolutionary models. 
Note that the depth to the 1.1887 fim Fe I line is overestimated in 
the model. 



Original SpeX + Best-Fit GAIA Model 
log(LB„/Lj = -3.37 ± 0.04 




Shifted SpeX + Best-Fit GAIA Model 
log(LB„/Lo) = -3.32 ± 0.05 




Fig. 15. — Bolomctric corrections for HD 114762B. The top 
panel shows the original SpeX/prism spectrum (dark grey) and 
the best-fitting model used for the bolomctric correction (light 
grey; [[M/H]/Teff /log g] = [-2.0/2500/1.0]). The bottom panel 
shows the SpeX spectral segments shifted to the published pho- 
tometry (dark grey) and the best-fitting model for the correction 
([+0.5/2600/3.0]). The models have been smoothed to the resolu- 
tion of the data. 



O'iRni = 47r(i^^Bo 




2ad 
d 



(8) 



The resulting bolometric flux and luminosity from the 
original spectrum are 9.11 ± 0.83 x 10^ W m~^ 



and log(LBoi/-^^o) = -3.37 ± 0.04, and for the shifted 
SpeX spectrum are 10.32 ± 0.95 x lO'^^ W m'^ and 
log(iBoi/^0) = -3.32 ± 0.04 (Figureim Tablep. 

We also investigated the effect of systematic errors 
due to the choice of model used in the bolometric 
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correction. We recalculated the bolometric luminosities 
of the original and shifted SpeX spectra using the 
second and third best-fitting synthetic spectra. The 
change in log(LBoi/i0) from choosing different model 
atmospheres was less than 0.01 dex for the original SpeX 
spectrum and was near 0.02 dex for the shifted spectrum. 
In addition, we computed the luminosities using the 
[-0.5/2600/5.5] GAIA model. This model was chosen 
because its parameters are the closest to those predicted 
by the evolutionary models (see § 3.3). The change in 
luminosiy was ^0.01 dex for the original spectrum and 
was ~ 0.02 dex for the shifted spectrum. These errors 
are uncorrelated for high signal-to-noise spectra, so we 
added the systematic errors to the statistical errors, in 
quadrature, to obtain our final uncertainty estimates in 
the luminosities: 0.04 and 0.05 dex for the original and 
shifted spectra, respectively. 



3.2.4. The Luminosities of Other Ultracool Subdwarfs 

Using the same method described in §3.2.3, we cal- 
culated the luminosities of four other ultracool sub- 
dwarfs with published parallaxes and whose spec- 
tra were available in the SpeX Prism Spectral L i- 
brary: LHS 377 rsdM7: IMonet et al.lfl99l. IGizisI 1199711. 



LSR 2036+5059 rsd M7.5: ILeoine et all I2003ir SSSPM 
1013-1356 (sdM9.5; [Scholz et alj [20Q4aD, and 2MASS 
iBure 



16264-3945 (sdL4; iBurgassen ,20041 ) . The spectra were 



flux calibrated using their 2MASS J-band magnitudes 
and the best-fitting GAIA models were used to de- 
termine bolometric corrections. Parallaxes for LSR 
2036+5059, SSSP M 1013-1356 , and 2MASS 1626+3945 

and the parallax for LHS 
The sp ectra were origi- 



are from Schilbach et al 
377 is from lMonet et al 




nallv pubhshed bv lBurgasseiT i 20041 ) and lBurgasser et all 
([^04). The best-fitting models used for the short- 
and long-wavelength corrections were [0.5/3200/5.5], [- 
1.0/2900/5.0]. [-1.5/2600/5.0], and [-2.5/2200/3.0] for 
LHS 377, LSR 2036+5059, SSSPM 1013-1356, and 
2MASS 1626+3945, respectively. We added in quadra- 
ture the resulting luminosity measurement errors and an 
estimated systematic error of 0.02 dex to account for the 
choice of atmospheric model used for the bolometric cor- 
rection. The results are summarized in Table [6l The 
luminosity of LHS 3 77 has previously been calculated by 
iLeggett et all (|2000D using an observationally-derived K- 
band bolometric correction; the luminosity we obtained 
(-3.08 ± 0.04 dex) is consistent with their value of -3.11 
dex. The luminosity of HD 114762B (-3.37 ± 0.04 dex) 
is similar to that of the sdM9.5 object SSSPM 1013-1356 
(-3.41 ± 0.08 dex), which is consistent with the spectral 
type of d/sdM9 ± 1 that we assigned to HD 114762B. 

In Figure \W\ we plot the bolometric luminosity 
as a function of opti c al spectral type for d warfs 
(from iDahn et all 120021 . iGolimowski et all 12004 and 



Gushing et al."2005') and subdwarfs (from this work and 
Burgasser et al. 2008b) . Objects that have since been 
resolved as close binaries have been removed from the 
samples. For a given optical spectral type, the luminosi- 
ties of ultracool subdwarfs appear to be slightly higher 
than those of dwarfs. A linear fit to the ultracool sub- 
dwarf luminosities from spectral types M7 to L7 yields 



-2.5 



-3.0 



o -3.5 



o -4.0 



-4.5 



-5.0 



Dahn et ol. (2002) Dwarfs * 
GolimowskI et ol. (2004) Dworfs ■ 
Gushing et ol. (2005) Dworfs O 
Subdwarfs • 




M6 



M8 



LO L2 L4 

Optical Spectral Type 



L6 



L8 



Fig. 16. — Bolometric luminosity as a function of spectral type 
for ultracool dwarfs (gray) and ultracool sub d warfs (red). The 
data f or the dwarfs are from IDahn et al.l II2002I ) , IGolimowski et ahl 
1I2OOI) . and ICushing et al.l II2005I) . Obiects that have since been 
discovered to be resolved close binaries have been removed from 
the samples. Th e subdw arf data arc from this work (Table |6ll and 
IBurgasser eTakl I l2008bl . for 2MASS 05325346+8246465 (sdL7)). 
The solid black line and the dashed black curves show the linear 
fit and ± 1 tr uncertainty (Equation 9). 



the following coefficients: 



log(L/LQ) = -2.258 - 0.111 x SpT, 
4.34 -0.369 



(9) 



Cov = 



-0.369 0.0357 



X 10" 



where SpT is the ultracool subdwarf numerical spectral 
type beginning at MO = and increasing by 1 for every 
spectral subclass (e.g., L5 = 15). Cov is the covariant 
matrix of the fit. 

There are several explanations for the apparent over- 
luminosity of ultracool subdwarfs in the luminosity- 
spectral type diagram. Unresolved binarity is one possi- 
bility, but it is unlikely that the few late-type subdwarfs 
sampled thus far all happen to be unresolved binaries. 
An alternative explanation is that the effective temper- 
atures and/or radii of ultracool subdwarfs are different 
than those of normal dwarfs in the same spectral sub- 
class. Ultracool subdwarfs appear to be ~ 0.2 dex above 
the luminosity-spectral type trend of ultracool dwarfs. If 
the higher luminosity were due to a larger radius alone 
then ultracool subdwarfs would need to have radii that 
were ~ 26% larger than those of dwarfs. If, on the 
other hand, the increased luminosity were caused by a 
difference in effective temperature, then ultracool sub- 
dwarfs would need effective temperatures that were ^ 
12% hotter than dwarfs from the same spectral sub- 
class. This translates into temperature differences of 
~ 150-250 K for objec ts throughout the late-L t o late - 
M spectral sequence. IBurgasser fc Kirkpatrickl ()2006f ) 
found additional evidence for different effective temper- 
ature scales for ultracool dwarfs and subdwarfs by com- 
paring their atmospheric model-derived effective temper- 
atures to their spectral types. They found that extreme 
M-type subdwarfs appear to have higher effective tem- 
peratures by several hundred Kelvin compared to normal 
M dwarfs, which they suggest might originate from dif- 
fering spectral classification schemes. Unequal effective 
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TABLE 5 

Luminosity and Radius of HD 114762B 







Ck (X 10-21) 


R {Re) 


R (Rq) 




Spectrum 


[[M/H]/T,ff/log g] 


{ = R^/<f} 


{ = dV^ } { = 


^LBol/(47r<7T4j) } 


log(LBol/io)" 


SpeX (1.00 - 2.35 Atm) 


[-2.0/2500/1.0] 


4.6 ± 0.4 


0.116 ± 0.006 


0.111 ± 0.010 


-3.37 ± 0.04 


Shifted SpeX {.J / H / K) 


[0.5/2600/3.0] 


4.4 ± 0.5 


0.113 ± 0.007 


0.110 ± 0.010 


-3.32 ± 0.05 



" Wc adopt a value of 3.86x10^'^ W m"^ for Lq. 



TABLE 6 

Luminosities and Bolometrk; Corrections of Ultracool Subdwarfs 



Object 


SpT 


Best-Fitting 
GAIA Model^' 


log(LBol/-f-0) 


BCj 
(mag) 


BCks 
(mag) 


Spectrum 
Ref. 


LHS 377 


sdM7 




[0.5/3200/5.5] 


-3.08 ± 0.04 


1.97 ± 0.04 


2.69 ± 0.04 


1 


LSR 2036+5059 


sdM7.5 




-1.0/2900/5.0 


-3.02 ± 0.05 


2.01 ± 0.04 


2.69 ± 0.04 


2 


SSSPM 1013-2356 


sdM9.5 




-1.5/2600/5.0 


-3.41 ± 0.08 


2.07 ± 0.05 


2.30 ± 0.09 


1 


2MASS 1626+3945 


sdL4 




-2.5/2200/3.0 


-3.71 ± 0.04 


2.21 ± 0.04 


2.18 ± 0.08 


1 



References. — f l) IBurgasserl I I2004I) . f2') IBurgasser et al.l I2004D 
^ For fits to the 0.80-2.35 /^m spectral region; model parameters are [[M/H]/T(.ff (K)/log g]. 



temperature scales for objects with the same spectral 
type but different luminosity class is not new: M dwarfs 
and M giants, for e xample, to exhibit this trait (see, e.g., 
Idi Benedettol[T99l . 

3.2.5. Bolometric Corrections for Ultracool Subdwarfs 

Using our computed bolometric fluxes, we can derive 
bolometric corrections for ultracool subdwarfs. The J- 
band bolometric correction is defined as BCj = AIbo\ - 
Mj = mBoi - J, where MboI is the absolute bolomet- 
ric magnitude, Afj is the J-band absolute magnitude, 
mBoi is the apparent bolometric magnitu de, and J is 
the J- band apparent magnitude. Following iBessell et al.l 
(|1998l Appendix D), we define MboI,© to be 4.74 and we 
assume Lbo\,q =3.86 x 10^^ W, in which case the J-band 
bolometric correction becomes BCj = — 2.5 log(i^Boi) — 
18.988 — J, where i^Boi is the bolometric flux in W. 

J- and ifg-band bolometric corrections for ultracool 
subdwarfs are plotted in Figure 1171 and are tabu l ated in 
Table m We use the values from lBurgasser et al.l ()2008bl ) 
for 2MASS 05325346+8246465 (sdL7) to derive bolomet- 
ric corrections for that object. A linear fit to the J-band 
bolometric corrections yield the following relation: 

BCj = 1.75 + 0.0328 x SpT, (10) 

„ _ f 4.95 -0.461 \ _,„_3 
Cov - y _o.461 0.0474 ) ^ ' 

where SpT is the numerical spectral type beginning at 
1VI7 = 7 and increasing by 1 for each spectral subclass 
through L7. The rms scatter about the fit is 0.007 mags. 
Similarly, a quadratic fit to the JiTs bolometric correction 
gives 

BCks = 5.05 - 0.461 X SpT + 0.0184 x S'pT^, (11) 

/ 292 -59.1 2.66 \ 
Cov = -59.1 12.0 -0.545 x 10"^ 
V 2.66 -0.545 0.0249 / 
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Fig. 17. — J- and Js's-band bolometric correction as a function 
of spectral type for ultracool subdwarfs. The gray dashed curve in 
the bottom panel shows the if - band bolometric co r rectio n (not /Cg ) 
for ultracool dwarfs derived bv lGolimowski et al.l l|2004l ). The dot- 
ted gray lines represent the typical rms scatter of ultracool dwarfs 
(0.13 mags) about the best-fitting polynomial relation. Ultracool 
subdwarf names and their individual corrections are tabulated in 
Table [6] The coefficients and covariance matrices of the fits are 
given in Equations 10 and 11. These relations are only applicable 
l^or mildly metal-poor "sd"-type subdwarfs. 

The rms scatter about the quadratic fit is 0.04 mags. We 
note that these relations are only applicable for moder- 
ately metal-poor ultracool subdwarfs and should not be 
applied for "usd"- or "esd"-type objects. 

3.2.6. The Radius of HD mi 6ZB 

The radius of HD 114762B can be calculated in two 
ways using the model atmospheres. One way is to use the 
scaling factor Cu which scales the emergent flux density 
of the best-fitting model to the fiux-calibrated spectra. 
The radius i? and the error aji can be computed using 
the following relations: 
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Fig. 18. — Predictions from the CB97 evolutionary models based on the bolometric luminosity of HD 114762B {log{L-Q^i/ Lq) = -3.37 
± 0.04, solid and dotted gray Unes) and the metallicity of HD 114762A {[Fe/H] = -0.71 ± 0.07). The [M/H] = -0.5 (red) and -1.0 (blue) 
metallicity traeks are shown for 1 Gyr (dotted lines) and 10 Gyr (solid lines). An interpolated —0.7 dex 10 Gyr metallieity track is shown 
in black. Mean values and errors (solid and dotted navy lines) for the effective temperature, radius, mass, and surface gravity are derived 
by performing a Monte Carlo simulation of luminosity and metallicity values. 



R^d^fcl, (12) 




where d and cr^ are the distance and uncertainty in the 
distance to the object, and Ck and a-^^ are the model 
scahng factor and its error. The radii computed using 
this method are 0.116 ± 0.006 Rq and 0.113 ± 0.007 Rq 
for the original and shifted SpeX spectra, respectively, 
where Rq = 6.960 x 10^° cm. 

The radius may also be computed using the bolometric 
luminosity and the effective temperature derived from 
the best-fitting atmospheric model: 

In this work we do not estimate the errors of the fun- 
damental parameters derived from fitting atmospheric 
models; instead we merely report the best few models 
as a triplet of [[M/H], T^a, log g\. To compute the un- 
certainty of the radius we estimate cTTeff to be ~ 100 K 



based on the narrow range of resulting temperatures from 
the best-fitting models. This gives radii of 0.111 ± 0.010 
Rq and 0.110 ± 0.010 Rq for the original and shifted 
SpeX spectra, respectively. These values are consistent 
with the radii derived using the model scaling factor. 

3.3. Comparing Physical Parameters From Atmospheric 
and Evolutionary Models 

Evolutionary models predict the change in stellar lu- 
minosity, effective temperature, and radius as a function 
of mass, age, and metallicity. When testing evolution- 
ary models, the interpretation of discrepancies between 
observations and model predictions must be made with 
care. Evolutionary models depend on atmospheric mod- 
els to account for the way that radiative flux escapes 
from a central luminosity source, so the choice of atmo- 
spheric models can influence the predictions of the evolu- 
tionary models. It may therefore be difficult to untangle 
inadequate atmospheric models from missing or incorrect 
physics in the evolutionary models. 

These models can be unraveled to some extent by 
studying their predictions in various diagrams. The 
fundamental physical parameters predicted by the evo- 
lutionary models (LboI, Toff, and R) are more depen- 
dent on t he input physics than on the atmospheric mod- 
els used (jChabrier et al.l l2000f ). However, converting 
L and T^e into observable quantities (magnitudes and 
colors) relies on the emergent spectral energy distribu- 
tions, so the predictive accuracy of color-magnitude di- 
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TABLE 7 

Results From CB97 Evolutionary Models for HD 
114762B 



Teff (K) R (Rq) M (Mq) logg 

2645 ± 42 0.100 ± 0.001 0.0879 ± 0.0007 5.381 ± 0.009 



Note. — Errors are derived in a Monte Carlo fashion using the 
bolometric luminosity from the original SpeX data (log(LBol/^o ) 
= -3.37 ± 0.04) and a motalhcity of -0.71 ± 0.07 dox. 

agrams reflects the accuracy of tlie atmosplieric mod- 
els more so than do L-Tcs HR diagrams. This depen- 
dency on correct atmospheric models is especiaUy strong 
for low-mass stars and brown dwarfs, for which the in- 
clusion of clouds in the models is required to even ap- 
proximately match the observed trends of late-M, L, 
and T dwarfs in the near-infrared color-absolute mag- 
nitude diagrams jA Uard c t alj 120011; iTsu ii & Na kaiimal 
[200l l&ilEows et ai.1.2006 : ISaumon fcMarlev 20081 ). In 
this study we use the fundame ntal parameters from by 
the iChabrier fc Baraffd (|T997l . CB 97) models a n d the 
color-magnitu de relations from th e iBaraffe et al.l () 19971 . 
BCAH97) and lBaraffe eFall (IT998I . BCAH98) models, all 
of which incorporate the same (PHOENIX/ A^eztGen) 
model atmospheres and are based on the same input 
physics. These models are widely used throughout the 
literature so it is crucial to test them at low masses and 
metallicities, a task that we perform here for the first 
time. 

Using the metallicity and the now well-measured lumi- 
nosity of HD 114762B, we compared the effective temper- 
atures, radii, and surface gravities from our atmospheric 
model fitting to the values predicted by the evolutionary 
models to test for consistency.^ Mutually consistent pre- 
dictions do not prove that the models are correct (both 
may predict the wrong values), but inconsistent results 
indicate that at least one set of models needs improve- 
ment. 

There are a large number of age estimates in the lit- 
erature for the primary star HD 114762A, the majority 
of which are based on evolutionary model isochrones or 
chromosphcric activity levels. Most estimates are greater 
than 10 Gyr (Table [T]). For the low mass evolutionary 
tracks that are relevant to HD 114762B, there is little 
difference between the 1 Gyr and 10 Gyr isochrones (Fig- 
ure [18]). We chose the 10 Gyr evolutionary tracks for our 
analysis. The properties of HD 114762B were derived in 
a Monte Carlo fashion by first randomly choosing a lu- 
minosity from a Gaussian distribution with mean values 
of -3.37 dex and a standard deviation of 0.04 dex. Simi- 
larly, a metallicity was chosen from a Gaussian distribu- 
tion with a mean value of -0.71 dex and a standard devi- 
ation of 0.07 dex. We created a new metallicity track at 
the chosen value by interpolating the evolutionary mod- 
els. Then using the luminosity and the new metallicity 
track, we determined a temperature, radius, mass, and 
surface gravity for that trial. This process was repeated 
10,000 times and the resulting distributions are summa- 
rized in Table El The CB97 models yielded the following 

^ While masses can be inferred from the atmospheric models 
using the radius and surface gravity, the results are not practical 
for comparison because the models are coarsely griddcd (Alog g = 
0.5). 



resuhs for log(LBoi/i0) = -3.37 ± 0.04 and [M/H] = - 
0.71 ± 0.07: Teff = 2645 ± 42 K, i? = 0.100 ± 0.001 Rq, 
M = 0.0879 ± 0.0007 Mq, and log g = 5.381 ± 0.009. 

The GAIA atmospheric models and the CB97 evolu- 
tionary models generally predicted consistent effective 
temperatures (between 2500 and 2700 K for the for- 
mer and 2645 K for the latter) and radii, but the in- 
ferred surface gravities were largely inconsistent. The 
exceptions were the fits to the OSIRIS data. Fitting 
the entire OSIRIS spectrum with the GAIA models re- 
sulted in a temperature of 2800 K (or 150 K above 
that of the CB97 evolutionary models) and a surface 
gravity of 5.5 dex (similar to 5.38 dex from the CB97 
evolutionary models). The fit to the J-band absorption 
lines y ielded T^s = 2600 K and log g = 5.0. iJao et all 
(l2008fl used the radii and masses of low-mass stars from 
iLopez-Mora lcs (2007') to construct a mass-gravity rela- 
tionship, which exhibited an increasing trend of surface 
gravity with lower mass. For a mass of ~ 0.09 Mq, we 
can therefore expect HD 114762B to have a surface grav- 
ity above ~ 5.0 dex. Most of the surface gravities from 
fits to the SpeX data were unrealistically low (log g < 
4.0) compared to the empirical trends and compared to 
the evolutionary models. Finally, we note that the ra- 
dius inferred from the atmospheric models (~ 0.11 Rq) 
is statistically slightly larger than the value predicted by 
the evolutionary models (~ 0.10 Rq). 

The color-magnitude diagram (CMD) provides another 
test of the atmospheric models. In Figure [19] we plot the 
position of HD 114762B with respect to the BCAH97 
and BCAH98 evolutionary models in the M^^g vs J - 
Ks CMD. The [M/H] = {0.0, -0.5} metalhcity tracks 
are from BCAH98 while the [M/H] = {-1.0, -1.3, -1.5, 
-2.0} tracks are from BCAH97. These evolutionary mod- 
els make use of an earlier version of the PHOENIX code 
(the NextGen models) than the GAIA grid that we stud- 
ied in this work. The location of HD 114762B in this 
diagram is comfortably between the -0.5 and -1.0 dex 
metallicity tracks, which is in precise agreement with the 
metallicity of -0.70 dex of the primary star. Previously 
wc showed that the synthetic spectra from the atmo- 
spheric models did a poor job of reproducing the correct 
physical parameters of HD 114762B from the best-fitting 
models. Nevertheless, the color-absolute magnitude dia- 
gram, which crudely measures the synthetic spectral en- 
ergy distributions of the models, appears to well-match 
the observations. This suggests that the shape of the 
models at lower metallicities is accurate, but that the de- 
tailed treatment of opacity sources in the near-infrared 
still needs improvement. 

To obtain a sense of the relative metallicities of other 
metal-poor low-mass objects, we also plot all currently 
known ultracool subdwarfs with parallax measurements 
in the Mk, vs J - Ks CMD in Fi gure [T9l Parallaxes 
are from iSchilbach et all (|2009[ ) and lMonet et all (|1992l . 
for LHS 377). Spectroscopically-classified mild subd- 
warfs fall in the expected mild-metallicity regions of the 
BCAH97/98 evolutionary tracks, and objects classified 
as subdwarfs are located between the [M/H] = -1.0 to 
-2.0 metallicity tracks. The M, L, and T spectral types 
appear to follow an increasing absolute magnitude trend 
with later spectral type, an expected feature for progres- 
sively later-type and lower-luminosity objects. The loca- 
tions of ultracool subdwarfs in the M^s vs J - Ks CMD 
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Fig. 19. — ^'^Kg vs J - Kg color-absolute magnitude diagram showing the positions of ultracool subdwarfs with parallaxes. The [M/H] 
= {0.0, -0.5} evolutionary tracks are from BCAH98 and the [M/H] = {-1.0, -1.3, -1.5, -2.0} tracks are from BCAH97 (10 Gyr models for 
both) . The evolutionary tra cks were converted fro m th e CIT photometric system to the 2M ASS system using the relations in ICarpenteii 
1I2OOII) . Parallaxes arc from lSchilbach et al.l I I2009I1 and liVlonei; et al.l (11993 '). See Figure 1 of lSchilbach et al.l (120091) for the identification 
of the other ultracool subdwarfs in this diagram. The position of HD 114762B in this diagram is shown as a yellow star. Its predicted 
metallicity from the color-magnitude diagram is consistent with the value of -0.70 dex for the primary star HD 114762A. 



therefore supports the current extension of the optical 
spectral classification scheme to lower metallicitics. Nev- 
ertheless, the development of a rigorous spectral classifi- 
cation scheme will be needed with the discovery of large 
numbers of ultracool subdwarfs. 

Based on the correct metallicity prediction for 
HD 114762B, we suspect that the accuracy of the 
BCAH97/98 evolutionary models is greater for ultracool 
subdwarfs than for their solar- metallicity counterparts, 
for which the inclusion of clouds in the atmospheric 
models is necessary to reproduce their observed J 
- Ks colors. In the low-mass metal-poor regime, 
continuum CIA H2 is the dominant opacity source in 
the near-infrared, so dust-free model atmospheres like 
NextGen probably do a reasonable job of reproducing 
the color-magnitude relations of ultracool subdwarfs. 
The discovery of additional benchmark companions over 
a large range of metallicitics will be necessary for this 
claim to be properly verified. 

3.4. Spectroscopic Parallaxes 
3.4.1. Spectroscopic Parallax to HD 114762B 

Since the factor Ck used to scale the model atmo- 
spheres to the flux calibrated spectra is equal to R^/cP, 
in principle accurate radius values could be used to de- 
rive the distance without knowing the parallax, assum- 
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Fig. 20. — Surface gravity as a function of effective tempera- 
ture for different values of radius and metal li city fr om the 10 Gyr 
evolutionary models of IChabrier fc Baraffel l|1997l ) . The models 
are interpolated over radius and metallicity. The gray shaded re- 
gion shows the range of T^ff from fitting synthetic spectra to the 
OSIRIS data (2800 ± 100 K). The red tracks show the metallicity 
range from the atmospheric model fitting (-0.5 ± 0.5 dex). The 
evolutionary models constrain the radius, which, combined with 
the model scaling factor Ck (= /d?), can be used to estimate 
the distance to HD 114762B without knowing the parallax or the 
bolometric flux and assuming the object is not a binary. 
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ing the object is not a binary. Evolutionary models can 
be used to constrain a radius if an age, effective tem- 
perature, and mctallicity are known. We have shown 
that, for HD 114762B, the most reliable technique for 
estimating the mctallicity, surface gravity, and effective 
temperature from fitting synthetic spectra is to use the 
the medium-resolution J-band region. Knowing a pos- 
teriori that this fitting technique is acceptable, we can 
use the resulting effective temperature and mctallicity to 
derive a radius from the evolutionary models and then 
a distance from Ck (equal to R^/d"^). This method is 
particularly well suited for observations of ultracool sub- 
dwarfs, which have ages that can be constrained from 
their kinematic properties. We refer to the distances es- 
timated using this method as "spectroscopic parallaxes." 

The best-fitting model to our OSIRIS spectrum is [- 
0.5/2800/5.5] (Table [3]), where, for this analysis, we as- 
sume errors of 0.5 dex for the mctallicity, 100 K for the 
effective temperature, and 0.5 dex for the surface grav- 
ity. For ultracool objects, the evolutionary models pro- 
vide tighter constraints on the surface gravity than the 
atmospheric models do, so the surface gravity is ignored 
for our estimate of the radius. The evolutionary models 
predict the radius to be between ~ 0.10 and 0.13 Rq for 
an effective temperature of 2800 ± 100 K, a mctallicity 
of -0.5 ± 0.5 dex, and an age of 10 Gyr, as shown in Fig- 
ure [20l We obtained a value of 3.82 ± 0.35 x lO^^i foj. 
Ck based on the best-fitting model to the fiux-calibrated 
OSIRIS data (§3.2.1). For a radius of 0.115 ± 0.015 i?©, 
the distance estimate is 42 ± 6 pc, where the errors have 
been propagated as in §3.2.5. The parallactic distance to 
HD 114762 A is 38.7 ± 1.1 pc. The spectroscopic paral- 
lax method using the OSIRIS data therefore predicts the 
distance to HD 114762B to within 10% of the parallactic 
distance. 

We note that this spectroscopic parallax technique is 
not strongly sensitive to effective temperatures cooler 
than ~ 3000 K. Although the distance is directly propor- 
tional to the radius, the tracks of constant radius tend to 
spread out at in log g-T^a space so that a small change 
in temperature produces a small change in radius in this 
low-temperature regime (Figure [201) • At hotter tempera- 
tures, however, a small change in temperature produces a 
large change in radius; the uncertainty in radius is there- 
fore impractically large for Toff ^ 3000 K. If we use a 
temperature of 2300 ± 100 K instead of 2800 K ± 100 
K for HD 114762B (to better match temperatures of late 
M dwarfs in the literature; see § 3.2.2), the resulting (ex- 
trapolated) radius is ~ 0.094 ± 0.009 (for a metaUicity 
of -0.5 ± 0.5 dex), and the distance estimate using the 
spectroscopic parallax technique for this cooler temper- 
ature is 34 ± 4 pc. This result is still consistent with the 
parallactic distance of HD 114762A at < 2 (t. 

3.4.2. Spectroscopic Parallaxes to Other Ultracool Subdwarfs 

We applied this same spectroscopic parallax technique 
to four ultracool subdwarfs with parallaxes previously 
discussed in §3.2.4 (LHS 377, LSR 2036+5059, SSSPM 
1013-1356, and 2MASS 1626+3945). As we have shown, 
the only reliable physical parameter from fitting atmo- 
spheric models to the SpeX spectra is the effective tem- 
perature. Nevertheless, it is still possible to obtain an 
estimate of the radius using the effective temperature by 
assuming metallicities less than solar but greater than 



TABLE 8 

Spectroscopic Parallaxes for Ultracool Subdwarfs 





Ck (x 10-21) 






7r 


Object 


{ = R'^/d'^ } 


(pc) 


(pc) 


Ref. 


LHS 377 


3.6 ± 0.1 


60 ± 30 


35 ± 1 


1 


LSR 2036+5059 


3.19 ± 0.09 


50 ± 10 


46 ± 3 


2 


SSSPM 1013-2356 


1.98 ± 0.06 


51 ± 8 


49 ± 5 


2 


2MASS 1626+3945 


4.1 ± 0.1 


32 ± 4 


34 ± 1 


2 



References. — d UMonet et al.llll992ll . (•2') ISchilbach et al.l I I2009D 
^ Spectroscopic parallax. See §3.4. 

1* Parallactic distance. Upper and lower limits on the distance errors 
are averaged for clarity. 



-2.0 dex (which encompasses the moderately mctal-poor 
"sd" class of ultracool subdwarfs), ages of ^10 Gyr. and 
errors of ± 100 K for the effective temperatures. Using 
the effective temperatures from the atmospheric model 
fits (Table[S]), we obtained radii estimates of 0.16 ± 0.07 
i?Q, 0.12 ± 0.03 i?0, 0.100 ± 0.015 Rq, and 0.090 ± 
0.010 Rq for LHS 377, LSR 2036+5059, SSSPM 1013- 
1356, and 2MASS 1626+3945, respectively. The result- 
ing distance estimates are listed in Table [8] along with 
the parallactic distances. The estimates from spectro- 
scopic parallaxes are surprisingly accurate, although less 
so at higher effective temperatures (> 3000 K) where the 
models spread out in surface gravity space and poorly 
constrain the radii values (Figure I20p . The precision in- 
creases at lower temperatures for the same reason^". 

Our spectroscopic parallax technique can be compared 
to the conventional photometric distance estimates which 
rely on empirical relations between spectral type and 
absolute magnitude. Polynomial fits to spectral type- 
absolute magnitude data for late-ty pe objects typically 
produce rms er rors of 0.3-0.4 mag (jTinnev et al.ll2003l : 
iLiu et alll2006[ ). which translate into distance errors of 
~ 15-20%. In comparison, our spectroscopic parallax 
method produces distances that are about twice as ac- 
curate (^ 10% agreement with published parallaxes), al- 
though this method has only been tested on a small sam- 
ple of subdwarfs. Extending this approach to a broader 
range of objects and wavelength ranges may be fruitful. 

4. SUMMARY AND CONCLUSION 

We have presented a near-infrared spectroscopic anal- 
ysis of HD 114762B, the first benchmark ultracool sub- 
dwarf companion. The mctallicity of HD 114762B is in- 
ferred from the primary star, which has a mean mctal- 
licity of [Fe/H] = -0.71 based on independent estimates 
from the literature. The spectral characteristics of HD 
114762B include suppressed H and A"-band continuum 
as well as deeper 1.4 and 1.9 /xm H2O bands compared 
to solar-metallicity M dwarfs. The SpeX 1.00-2.35 fim 

1" We note that ILeggett et al.l II2000I ) obtained an effective tem- 
perature of 2900 ± 100 K for LHS 377 using TaB-{V-I)-[M/U] 
relations from evolutionary models. A visual inspection of the best- 
fitting GAIA model shows poor agreement between the LHS 377 
spectrum and the model, so our value of 3200 K from the best- 
fitting GAIA model may be overestimated. A reduced T^g for 
LHS 377 would lower the radius obtained from the evolutionary 
models, which would place the estimated distance closer to the ac- 
tual value. This may suggest that GAIA model fits to earlier-type 
ultracool subdwarfs are not reliable, although we have shown that 
consistent values were obtained for objects near d/sdM8-d/sdM9 
types (§3.2.2). 
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spectrum does not, however, exhibit the extreme CIA H2 
which is characteristic of the "sd" class of ultracool sub- 
dwarfs. In this respect, HD 114762B is more similar to 
mild subdwarfs, which are believed to have an intermedi- 
ate metallicity between dwarfs and subdwarfs. We assign 
it a spectral type of d/sdM9 ± 1 based on its resemblance 
to 2MASS 0041-^35 (d/sdM9) from 1.00 to 2.35 fim. A 
spectral type later than sdM7.5 can also be inferred from 
the strength of the K I lines and the strength of the 1.4 
fim H2O band compared to M-type subdwarf spectra. 
The measured luminosity places HD 114762B just above 
the hydrogen-burning minimum mass based on the evo- 
lutionary models of CB97. 

We use this unique benchmark object to test low-mass 
metal-poor atmospheric and evolutionary models. The 
best-fitting PHOENIX/GA/^ synthetic spectrum to our 
OSIRIS J-band data has physical parameters that are 
consistent with both the metallicity of the primary star 
and the high surface gravity of late-type objects. Fits to 
the SpeX 1.00-2.35 ^m spectrum, however, yield physical 
parameters that are inconsistent with these values. We 
also fit the models to low-resolution spectra of other late- 
type mild ultracool subdwarfs with similar spectral types 
to HD 114762B. The best-fitting models' surface gravity 
and metallicity are mutually inconsistent among the four 
objects and were likewise inconsistent for fits to differ- 
ent regions of the same spectrum. Fits to HD 114762B, 
however, yield consistent effective temperatures between 
2500-2800 K. We conclude that the metallicities and sur- 
face gravities derived from fitting atmospheric models 
to low-resolution near-infrared spectra of ultracool sub- 
dwarfs are not trustworthy, but that the effective tem- 
peratures are probably more reliable. 

Problems with atmospheric models at low masses and 
low metallicities have already been noted in the litera- 
ture; only moderately good matches have resulted from 
fits to the red-optical spectral regions of ultracool sub- 
dwarfs, even when applied over rela.tively short wave- 
length rang es ("Sch weitzer et al.|[l999l : iLepine et al.ll2004l : 
iBurgasser c t al. 2007f). These same studies held the 
model surface gravity fixed at log g = 5.0 or 5.5 while 
only allowing the temperature and metallicity to vary. 
The models should, however, be able to predict all three 
parameters with no a priori assumption about any par- 
ticular one. In addition, the models should be able to 
predict roughly the same parameters for objects of the 
same spectral type as well as for different spectral regions 
of the same object. We have shown that neither appears 
to be true for the ultracool subdwarfs we studied, mak- 
ing the discovery of benchmark metallicity calibrators 
like HD 114762B ever more important. 

Our results suggest that the best method to determine 
the metallicity of field ultracool subdwarfs may be to 
use the medium-resolution J-band region. The 1.18-1.35 
fim range encompasses metal-sensitive atomic absorp- 
tion lines as well as the metal-sensitive 1.4 fim steam 
band. An interesting line that may be important for 
very metal-poor ultracool subdwarfs is the Fe I 1.887 
/im feature, which the GAIA models predict to be es- 
sentially independent of temperature and gravity, but 
which is strongly dependent on metallicity for [M/H] < 
-1.0. We found that fitting the individual absorption 
lines rather than the entire 1.18-1.35 /im spectrum of 
HD 114762B produced temperatures and gravities that 



were more consistent with the values derived from the 
CB97 evolutionary models. Nevertheless, fitting the en- 
tire J-band spectrum and the individual absorption fea- 
tures yielded metallicities consistent with the primary 
star and surface gravities consistent with empirical val- 
ues of late-type objects. 

The location of HD 114762B in the BCAH97/98 M^^ 
vs J - JsTs color-absolute magnitude diagram is consistent 
with the metallicity of the primary star. The metallic- 
ity and surface gravity inferred from GAIA fits to low- 
resolution near-infrared spectra are not trustworthy, but 
the results from the color-magnitude diagram suggest 
that the overall shape of the synthetic spectra are in 
fact reliable, at least for the mass 0.088 Mq) and 
metallicity ([Fe/H] = -0.71) of HD 114762B. 

Additionally, we calculated the bolometric luminosities 
and bolometric corrections for four known ultracool sub- 
dwarfs (LHS 377, LSR 2036+5059, SSSPM 1013-1356, 
and 2MASS 1626-f3945) with parallaxes and whose near- 
infrared spectra were available in the SpeX Prism Spec- 
tral Library. Ultracool subdwarfs have higher luminosi- 
ties and smaller bolometric corrections than do ultracool 
dwarfs in the same spectral subclass. While unresolved 
binarity is a possibility, a hotter effective temperature 
scale and/or a larger radius is a more likely explanation 
for the overluminosity. 

Finally, we have also developed a technique to estimate 
the distances to ultracool subdwarfs based on the avail- 
able models and a single near-infrared spectrum. Using 
the effective temperature from the best-fitting synthetic 
spectrum and an assumed age (~ 10 Gyr for objects with 
thick disk or halo kinematics), a radius can be inferred 
from evolutionary models. The radius can then be used 
along with the model scaling factor Ck (equal to /cP) 
to derive a distance estimate to the object. We applied 
this technique to five ultracool subdwarfs with trigono- 
metric distances and obtained distance estimates accu- 
rate to < 10% of the parallactic distances, or about 2 
times better than estimates based on photometry alone. 
This technique is particularly useful for ultracool subd- 
warfs, whose ages can be constrained by their kinematics. 

It will be worthwhile to study HD 114762B at other 
wavelengths, including obtaining thermal infrared pho- 
tometry (A > 3 ^m) to better constrain the long- 
wavelength end of the model-fitting, and, if possible, 
optical spectra, which would allow for an independent 
spectral classification of this object. The dominant error 
in our luminosity determination originates from the er- 
rors in the photometry used to fiux calibrate the spectra. 
Better near-infrared photometry would therefore enable 
more a precise luminosity estimate and, by extension, 
more precise predictions from the evolutionary models. 
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National Science Foundation. MCL and BPB acknowl- 
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Fellowship, NSF grant AST-0507833. The authors wish 
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